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Abstract. We study a general class of nonlinear mean field Fokker-Planck equations in relation with an 
effective generalized thermodynamical (E.G.T.) formalism. We show that these equations describe several 
physical systems such as: chemotaxis of bacterial populations, Bose-Einstein condensation in the canoni- 
cal ensemble, porous media, generalized Cahn-Hilliard equations, Kuramoto model, BMF model, Burgers 
equation, Smoluchowski-Poisson system for self-gravitating Brownian particles, Debye-Hiickel theory of 
electrolytes, two-dimensional turbulence... In particular, we show that nonlinear mean field Fokker-Planck 
equations can provide generalized Keller-Segel models for the chemotaxis of biological populations. As an 
example, we introduce a new model of chemotaxis incorporating both effects of anomalous diffusion and 
exclusion principle (volume filling). Therefore, the notion of generalized thermodynamics can have appli- 
cations for concrete physical systems. We also consider nonlinear mean field Fokker-Planck equations in 
phase space and show the passage from the generalized Kramers equation to the generalized Smoluchowski 
equation in a strong friction limit. Our formalism is simple and illustrated by several explicit examples 
corresponding to Boltzmann, Tsallis, Fermi-Dirac and Bose-Einstein entropies among others. 
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1 Introduction 

The study of Fokker-Planck equations (Fokkcr (1914) p], 
Planck (1917) [2]) is of considerable interest in physics 
since the pioneering work of Einstein (1905) [3] on the 
Brownian motion. The simplest Fokker-Planck equation 
is the Smoluchowski (1915) [4] equation. This is a drift- 
diffusion equation describing the diffusion of particles in 
physical space in the presence of an external potential 
(e.g. the sedimentation of colloids in a gravitational field). 
A more general Fokker-Planck equation is the Kramers 
(1940) [5j equation, previously introduced by Klein (1921) 
[BJ, which takes into account inertial effects and describes 
the diffusion of particles in phase space when they ex- 
perience a friction force. The Smoluchowski equation is 
recovered from the Kramers equation in a strong friction 
limit where inertial effects are negligible. In this sense, 
the Smoluchowski equation describes an overdamped evo- 
lution. These Fokker-Planck equations [7j are consistent 
with usual thermodynamics in the canonical ensemble. 
They monotonically decrease the Boltzmann free energy 
(iJ-theorem) and relax towards the equilibrium Boltzmann 
distribution. 

Recently, several researchers have tried to extend the 
usual concepts of thermodynamics and kinetic theory in 
order to describe complex systems that arc characterized 



by non-Boltzmannian distributions. In that respect, some 
generalized forms of entropic functionals Q have been in- 
troduced. One of the most popular "generalized entropy" 
is the Tsallis (1988) [3] entropy, but other entropies have 
been presented by Abe (1997) [TU], Borges & Roditi (1998) 
HQ, Kaniadakis (2001) pj], Naudts (2004) [13], and Ka- 
niadakis et al. (2005) [T3]. It was later realized that these 
entropic functionals arc special cases of the one-parameter 
family of entropies introduced earlier by Harvda & Char- 
vat (1967) [T5] or of the two-parameters family of en- 
tropies introduced by Mittal (1975) [16] and Sharma & 
Taneja (1975) [17] ■ Other famous forms of entropies have 
been presented by Reyni (1970) [18] and Sharma & Mittal 
(1975) [19]. We refer to Kaniadakis & Lissia [20] for a very 
interesting discussion of these historical aspects, starting 
from the early works of Euler in 1779. 



We recall that the Boltzmann entropy can be obtained 
from a combinatorial analysis assuming that all the accessi- 
ble microstates are equiprobable. This is the basic postulate 
of statistical mechanics. Non-standard entropies can be rele- 
vant for complex systems where this postulate breaks down, 
i.e. when the accessible microstates are not equiprobable. This 
happens when the system prefers some regions of phase space 
better than others or when the particles are subjected to ex- 
clusion/inclusion principles or fine-grained constraints [8]. 
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Following these developments, some researchers have 
tried to develop out-of-equilibrium theories associated to 
a generalized thermodynamical framework. In particular, 
it has been first shown by Plastino & Plastino (1995) [21] 
that the Tsallis q-distributions are the steady states of 
a nonlinear Fokker-Planck equation taking into account 
anomalous diffusion. This type of equations had been pre- 
viously considered by mathematicians to describe porous 
media [22] . The seminal work of Plastino & Plastino [2T] 
has been further developed by Tsallis & Bukman (1996) 
[23], Stariolo (1997) [2J, Borland (1998) [25] and Nobre 
et al. (2004) [26] among others. On the other hand, Ka- 
niadakis & Quarati (1994) [27] have introduced nonlin- 
ear Fokker-Planck equations whose steady states are the 
Fermi-Dirac and Bose-Einstein statistics. These kinetic 
equations take into account an exclusion (fermions) or in- 
clusion (bosons) principle leading to quantum-like statis- 
tics at equilibrium. The case of intermediate statistics, 
interpolating between fermions and bosons, has also been 
considered in [27j . Recently, the bosonic Kramers equa- 
tion has been studied in [30] and was shown to reproduce 
the phenomenology of the Bose-Einstein condensation in 
the canonical ensemble. 

The above-mentioned nonlinear Fokker-Planck (NFP) 
equations are associated with special forms of cntropic 
functionals (Tsallis, Fermi-Dirac, Bose-Einstein). More re- 
cently, Martinez et al. (1998) [31], Kaniadakis (2001) [12], 
Frank (2002) [32] and Chavanis (2003) [33] have studied 
generalized forms of NFP equations associated with an al- 
most arbitrary cntropic functional. They can be viewed as 
generalized Kramers and Smoluchowski equations where 
the coefficients of diffusion, friction and drift explicitly 
depend on the local density of particles. Physically, this 
can take into account microscopic constraints (exclusion 
volume effects, stcric hindrance, non-extensive effects...) 
that modify the dynamics of the particles at small scales 
and lead to non-standard equilibrium distributions 0. Mar- 
tinez et al. [31] determined the NFP equation in order 
to recover, as a steady state, the equilibrium state pro- 
duced by minimizing a generalized form of free energy at 
fixed mass. Kaniadakis Q~2] obtained the NFP equation 

2 A generalized Fokker-Planck equation leading to the Fermi- 
Dirac statistics has also been introduced by Chavanis et al. 
(1996) 28 in the context of the violent relaxation of collision- 
less stellar systems described by the Vlasov equation. This is 
based on the Lynden-BelPs form of entropy (1968) [2S] which 
becomes similar to the Fermi-Dirac entropy in the two-levels 
approximation of the theory. 

3 Generalized Kramers and Smoluchowski equations describe 
dissipative systems where the temperature is fixed instead of 
the energy. They are therefore associated with the canonical 
ensemble. The appropriate thermodynamical potential is the 
free energy F = E — TS which decreases monotonically with 
time at fixed mass. Generalized Boltzmann and Landau equa- 
tions describing conservative systems where the energy is fixed 
have been introduced by Kaniadakis (2001) [12] and Chavanis 
(2004) [34] . They are associated with the microcanonical en- 
semble. The proper thermodynamical potential is the entropy 
S which increases monotonically with time at fixed mass and 
energy. 



from the Master equation (see also Curado & Nobre [35] ) 
by allowing the transition probabilities to depend on the 
concentration of particles in the initial and arrival states. 
Frank [32] derived the NFP equation from a generalized 
free energy functional by using the linear thermodynamics 
of Onsager. Chavanis [33] obtained the NFP equation by 
using a form of Maximum Entropy Production Principle 
(MEPP). This corresponds to a variational version of the 
linear thermodynamics of Onsager. We refer to the book 
of Frank [36] for a first survey on nonlinear Fokker-Planck 
equations. 

Another topic of active research in statistical physics 
concerns the dynamics and thermodynamics of systems 
with long-range interactions |37] . Several Hamiltonian sys- 
tems with long-range interactions have been studied in 
the microcanonical ensemble such as self-gravitating sys- 
tems, two-dimensional vortices, the HMF model, the free- 
electron laser etc... For such systems the mean field ap- 
proximation provides a very good description of the sys- 
tem and becomes exact in a proper thermodynamic limit 
N — > +oo. The microcanonical ensemble is the correct de- 
scription of isolated systems evolving at fixed energy. On 
the other hand, some authors have introduced a canon- 
ical version of these models so as to treat systems that 
are dissipative. This leads to the notion of Brownian sys- 
tems with long-range interactions. Their dynamics is de- 
scribed by mean field Fokker-Planck equations where the 
temperature is fixed (instead of the energy). These mean 
field Fokker-Planck equations were introduced early by 
Kuramoto (1984) [38] to describe the synchronization of 
globally coupled nonlinear oscillators and more recently 
by Marzel & Aslangul (2001) [39], Chavanis (2006) [40] 
and Frank |36] in a more general context. Some specific 
studies have been made for self-gravitating Brownian par- 
ticles |41|42j and for the BMF model [43] which is the 
canonical version of the HMF model [44] . 

In view of the importance of these two topics: general- 
ized thermodynamics and long-range interactions, we have 
introduced in [33] a class of nonlinear mean field Fokker- 
Planck equations (see Eq. (81) of [33]) that incorporate 
both a generalized free energy functional and a long-range 
potential of interaction. As an illustration, we studied with 
C. Sire a model of self-gravitating Langevin particles [?5] 
(see also [46]) that combines self-gravity (long-range in- 
teractions) and anomalous diffusion (generalized thermo- 
dynamics) related to the Tsallis entropy and to the poly- 
tropic equation of state. In other words, this model couples 
the NFP equation introduced by Plastino & Plastino [21] 
to the gravitational Poisson equation. Unfortunately, this 
model of self-gravitating Langevin particles has no clear 
application in astrophysics (because self-gravitating sys- 
tems are generally not overdamped and not dissipative) 
and was introduced essentially as an interesting dynami- 
cal model with rich mathematical properties. However, it 
was noted in [45] that this model could have applications 
in unexpected area, like in the chemotaxis of bacterial 
populations... 

The name chemotaxis refers to the motion of organ- 
isms induced by chemical signals |47j . In some cases, the 
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biological organisms (bacteria, amoebae, endothelial cells, 
ants...) secrete a substance (pheromone, smell, food, ...) 
that has an attractive effect on the organisms themselves. 
Therefore, in addition to their diffusive motion, they move 
preferentially along the gradient of concentration of the 
chemical they secrete (chemotactic flux). When attrac- 
tion prevails over diffusion, the chcmotaxis can trigger 
a self-accelerating process until a point at which aggre- 
gation takes place. This is the case for the slime mold 
Dictyostelium discoideum and for the bacteria Escherichia 
coli. This is referred to as chemotactic collapse. A model 
of slime mold aggregation has been introduced by Pat- 
lak (1953) [48] and Keller & Segel (1971) [49] in the form 
of two coupled differential equations. The first equation 
is a drift-diffusion equation describing the evolution of 
the concentration of bacteria and the second equation is 
a diffusion equation with terms of source and degrada- 
tion describing the evolution of the concentration of the 
chemical. In the simplest model, the diffusion coefficient 
D and the mobility x of the bacteria are constant. This 
forms the standard Keller-Segel model. However, the orig- 
inal Keller-Segel model allows these coefficients to depend 
on the concentration of the bacteria and of the chemical. 
The case where these coefficients depend on the concen- 
tration of the chemical c(r, t), but not on the concentra- 
tion of bacteria p(r, t) , has been considered by Othmer 
& Stevens [50] . This leads to ordinary mean field Fokker- 
Planck equations (with respect to p(r,t)) with space and 
time dependent coefficients. On the other hand, if we as- 
sume that the diffusion coefficient and the mobility of the 
bacteria depend on their concentration p(r, t) , but not on 
the concentration c(r, t) of the secreted chemical, the orig- 
inal Keller-Segel model takes the form of a generalized 
mean field Fokker-Planck equation. Therefore, the Keller- 
Segel model represents a fundamental example of nonlin- 
ear mean field Fokker-Planck equation with physical ap- 
plications in biology, thereby justifying a notion of gener- 
alized thermodynamics. The analogy between the Keller- 
Segel model and NFP equations was first pointed out in 
[33 and further developed in subsequent papers (see, e.g., 
EI]). This analogy makes possible to interprete results of 
chcmotaxis in terms of nonlinear Fokker-Planck equations 
and generalized thermodynamics, which has not been con- 
sidered so far by applied mathematicians working on this 
domain |52] . It thus allows to make a bridge between two 
different communities. In this analogy, the model of self- 
gravitating Langevin particles introduced by Chavanis & 
Sire [45] also provides a generalized Kcllcr-Scgcl model of 
chemotaxis taking into account anomalous diffusion (the 
application of this model to chemotaxis has been empha- 
sized in [53]). More generally, we can use the numerous 
results accumulated in the context of generalized thermo- 
dynamics to propose new forms of generalized Keller-Segel 
models with potential applications in biology. 

The aim of this paper is to develop a simple and rich 
formalism that allows to deal with nonlinear mean field 
Fokker-Planck equations. We shall illustrate this formal- 
ism on several examples and show the inter-connections 
between different topics. The paper is organized as follows. 



In Sec. [21 we consider NFP equations in physical space. 
This corresponds to overdamped models where inertia! 
effects are neglected. We review and complete the basic 
properties of these equations in relation with an effective 
generalized thcrmodynamical (E.G.T.) formalism. In Sec. 
12. H we show that they can be obtained from generalized 
Langevin equations. In Sec. 12. 5( we show that they admit 
an .ff-theorem (in the canonical ensemble) for a general- 
ized free energy. We stress that the Legendre structure of 
the free energy and the Einstein relation are preserved in 
this generalized thcrmodynamical framework. In Sec. 12 . Gl 
we determine the steady states of these nonlinear mean 
field Fokker-Planck equations and show that they are so- 
lutions of an integrodiffcrcntial equation. In Sec. 12.71 (and 
in Appendix IE]), we show that a steady state of a nonlinear 
mean field Fokker-Planck equation is linearly dynamically 
stable if and only if (iff) it is a (local) minimum of the 
free energy at fixed mass. In Sec. 12. 9[ we show that a 
NFP equation in physical space with a constant mobility 
and a density-dependent diffusion coefficient can be writ- 
ten in the form of a generalized Smoluchowski equation 
incorporating a barotropic equation of state. In Sees. 12.101 
and 12.111 we show the corrcspondancc between the phc- 
nomenological derivations of the NFP equations given by 
Frank [32] and Chavanis [33] and the kinetic derivation 
given by Kaniadakis [12j . In Sec. El we present several ex- 
plicit examples of NFP equations and mention their po- 
tential applications to the problem of chemotaxis (see Sec. 
12. 4p . In particular, we introduce a new model of chemo- 
taxis that incorporates both effects of anomalous diffusion 
and exclusion principle (volume filling). The correspond- 
ing generalized entropy is expressed in the form of inte- 
grals that can be explicited in particular cases. In Sec.[U 
we consider NFP equations in phase space taking into ac- 
count incrtial effects. In Sees. I4.HI4.7[ we list their main 
properties. In Sec. 14.81 we consider the strong friction 
limit and derive the generalized Smoluchowski equation 
from the generalized Kramers equation. We use a method 
of moments that is simpler than the Chapman-Enskog 
expansion presented in [54] . In Sec. El we consider ex- 
plicit examples corresponding to the Boltzmann, Tsallis 
and Fermi-Dirac entropies. The Appendices contain im- 
portant results that complete the basic properties of the 
NFP equations discussed in the text. In Appendix |A1 we 
show that a generalized isotropic BGK operator has prop- 
erties similar to those possessed by a nonlinear Kramers 
operator. In Appendix [Bl we establish a simple relation 
showing the equivalence between linear dynamical stabil- 
ity (exponential damping of the perturbation) and gener- 
alized thcrmodynamical stability (minimum of free energy 
at fixed mass). In Appendix [C] we study the stability of a 
spatially homogeneous solution of the nonlinear mean field 
Fokker-Planck equation and evidence a critical point. In 
Appendices [Pi and lEl we show that a distribution function 
/ in phase space is a minimum of the free energy F[f] at 
fixed mass iff the corresponding distribution p in physical 
space is a minimum of the corresponding free energy F[p] 
at fixed mass. This implies that a distribution function / 
in phase space is a stable steady state of the generalized 
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Kramers equation iff the corresponding distribution p in 
physical space is a stable steady state of the correspond- 
ing generalized Smoluchowski equation. In Appendix [Fj 
we extend to d dimensions the Eddington formula that al- 
lows to obtain the distribution function / = /(e) in phase 
space from the knowledge of the barotropic equation of 
state p = p(p) in physical space. In Appendix [Gj we de- 
rive the if-theorems associated with the NFP equations 
and in Appendix [H] we derive the polytropic equation of 
state associated with the Tsallis statistics. 



2 Nonlinear mean field Fokker-Planck 
equations in physical space 

We first describe nonlinear mean field Fokker-Planck equa- 
tions in physical space where the inertia of the particles 
is neglected. They can be viewed as models describing an 
overdamped dynamics. 

2.1 Generalized Langevin equations 

We consider a system of ./V particles whose individual dy- 
namics is described by the stochastic Ito-Langcvin equa- 
tions 

dr 

= -x(Pi)V$i + v/25fo)Ri(*). (!) 



where Rj(i) is a white noise satisfying (Rj(t)) = and 
(Ri, a (t)Rj,p(lf)) = 6ijS a!0 5(t - t') where i = 1,...,N la- 
bel the particles and a = 1, ...,d label the coordinates of 
space. We have noted pi = p(ri(t), t) and = $(ri(t), t). 
In ordinary models, the mobility \ an d the diffusion co- 
efficient D are constant. In that case, the statistical equi- 
librium state is the Boltzmann distribution p ~ e~*/ T 
where the temperature T = 1//3 is given by the Einstein 
relation T = D/x- In the present study, we shall consider 
more general situations and allow the mobility x{p) an d 
the diffusion coefficient D(p) to depend on the local con- 
centration of particles p(v,t) = Q^j=i ^( r — r iW))- This 
is an heuristic approach to take into account microscopic 
constraints that affect the dynamics of particles at small 
scales and lead to non-Boltzmannian equilibrium distri- 
butions. Indeed, it is not surprising that the mobility or 
the diffusive properties of a particle depend on its envi- 
ronment. For example, in a dense medium, its motion can 
be hampered by the presence of the other particles so that 
its mobility is reduced. 

On the other hand, in ordinary models, the particles 
move in a fixed external potential <P ext (r). In the present 
study, we want to be more general and take into account 
the possibility that the potential ^(r, t) is created self- 
consistently by the particles themselves. In this paper, we 
shall neglect statistical correlations and use a mean field 
description (for more general models taking into account 
statistical correlations sec, e.g., |4QI55j ). Therefore, we as- 
sume that the potential is given by a relation of the form 



where u(\r — r'|) is a binary potential of interaction and 
p(r, t) is the smooth distribution of particles. In general, 
the mean field approximation gives a very good descrip- 
tion of systems with weak long-range binary interactions 
and it becomes exact in a proper thermodynamic limit 
N — * +oo [40]. In Eq. ([2]), the potential is expressed as a 
convolution product: <P = u * p. Of course, the potential 
can be due to the combined effect of the self-interaction 
and an external field, in which case = <P e xt +u* p. We 
shall also consider the case where it is determined by an 
equation of the form 



= A$ - k 2 $ - Xp, 

ot 



(3) 



where e and A are positive constants. For e = 0, Eq. ^ 
becomes the screened Poisson equation 

Z\<Z> - fc 2 <Z> = Xp. (4) 

Therefore, we can identify k^ 1 as the screening length. If 
we assume furthermore that k — 0, we get the Poisson 
equation 



A$ = Xp. 



(5) 



These last two equations can be put in the form of Eq. . 
Note also that, in the stationary state, Eq. ([3]) reduces to 
Eq. Q. 



2.2 Drift-diffusion equations 

For the stochastic process ([1]), the evolution of the smooth 
density of particles p(r, t) is governed by the nonlinear 
mean field Fokker-Planck equation [33I40| : 



dp 
dt 



V ■{V(D(p)p)+ X (p)pV<P}, 



(G) 



coupled to Eq. (f2|) or (J3j) . Let us introduce the notations 



Dh(p) = —(pD(p)), 
dp 



X9(j>) = PX(p), 



(T) 



where D and x are positive constants and h(p) and g(p) 
are positive functions. These notations are chosen such 
that the usual stochastic equations with constant diffusion 
D{p) — D and constant mobility x{p) = X are recovered 
for h(p) = 1 and g(p) = p. With these notations, the 
nonlinear Fokker-Planck equation ([6]) can be rewritten 



= V • (Dh(p)Vp + X9(P)V&) 
It can be put in the conservative form 

where 



(8) 



(9) 



«P(r,*)= / p(r',t)u(\r-r'\)dr', 



(2) 



J = -(Dh(p)Vp + X9(p)V$), (10) 

is a diffusion current. This structure guarantees the con- 
servation of mass M = J pdr provided that the normal 
component of the current at the boundary vanishes. 
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2.3 Relation to some known models 

Known models can be recovered as particular cases of 
Eq. (J6j) . When <P ex t is an external potential and when 
D{p) = D and x{p) — X arc constant, we recover the 
Smoluchowski equation d t p = V ■ {DV ' p + xP^^ext) de- 
scribing the Brownian motion of colloidal suspensions in 
a fixed gravitational field [4]. When pD(p) = Kp 1 and 
^ext = 0, we recover the porous medium equation dtp = 
KAp 1 [22], and when <P ex t ^ 0, we recover the nonlinear 
Fokker-Planck equation d t p = V • (i^Vp 7 + xP^^ext) in- 
troduced by Plastino & Plastino [21] in connection to the 
Tsallis statistics [9]. When D(p) = D and x(p) = X are 
constant and u = —^cos(9 — 8') we obtain the Brown- 
ian Mean Field (BMF) model which is the canonical 
(fixed T) version of the microcanonical (fixed E) Hamil- 
tonian Mean Field (HMF) model [44]. This is closely re- 
lated to the Kuramoto model [3BJ which describes the 
synchronization of globally coupled nonlinear oscillators. 
When D(p) = D and x(p) = X are constant and <P is the 
gravitational potential satisfying A<P = AirGp, Eqs. ([6j- 
((5]) become the Smoluchowski-Poisson system describing 
a model of self-gravitating Brownian particles studied by 
Chavanis & Sire [41142156157158159] . When D(p) = D and 
x(p) = X are constant and & is the self-consistent electro- 
static potential, we recover the equations introduced by 
Dcbyc & Huckel [60] in their model of electrolytes. Models 
of the form ([6])- ([4]) including a time dependent tempera- 
ture (3(t) assuring the conservation of energy have been 
introduced by Robert & Sommeria [61] and Chavanis [62] 
to describe the violent relaxation of two-dimensional vor- 
tices in geophysical and astrophysical flows. In the Quasi 
Geostrophic (QG) model, /c _1 represents the Rossby length 
Similar equations have been proposed in |28] to model the 
violent relaxation of stellar systems. On the other hand, 
for short range interactions, we can expand the poten- 
tial in the form <£(r,i) = ap(r,i) + |Z\/o(r,t) and Eq. 
([6j leads to a generalized form of the Cahn-Hilliard equa- 
tion (see [54151155] for details). As a particular case, for 
D = and u — a5 (leading to $(r,t) = ap(r,t)), we get 
d t p = %aV(g(p)Vp) and for g(p) = p. we get the porous 
medium equation dtp = \x&Ap 2 . Finally, in d = 1, as- 
suming D(p) = D, x(p) = X an d u' = -^5, we get 
the Burgers equation d t p + pd x p = Dd xx p [53]. Therefore, 
the class of nonlinear mean field Fokker-Planck equations 
([6j-([2]) introduced in |33] can find physical applications in 
different areas [5]- 



2.4 Generalized Keller-Segel model of chemotaxis 

In addition to the previous examples, nonlinear mean field 
Fokker-Planck equations can find important applications 
in the context of chemotaxis [47] . The original Keller-Segel 
model [49] describing the chemotaxis of bacterial popula- 
tions consists in two coupled differential equations 

% = V • (D 2 V P ) - V ■ (AVc) , (11) 



e— = D c Ac-k(c)c + f(c)p, (12) 

that govern the evolution of the density of bacteria p(r, t) 
and the evolution of the secreted chemical c(r, t) . The bac- 
teria diffuse with a diffusion coefficient D% and they also 
move in a direction of a positive gradient of the chemi- 
cal (chemotactic drift). The coefficient D\ is a measure of 
the strength of the influence of the chemical gradient on 
the flow of bacteria. On the other hand, the chemical is 
produced by the bacteria with a rate /(c) and is degraded 
with a rate fc(c). It also diffuses with a diffusion coefficient 
D c . In the general Keller-Segel model, D\ = Di(p, c) and 
Z?2 = D2 {p, c) can both depend on the concentration of 
the bacteria and of the chemical. This takes into account 
microscopic constraints, like close-packing effects, that can 
hinder the movement of bacteria. If we assume that the 
quantities only depend on the concentration of bacteria [f] 
and write D 2 = Dh(p), D\ = x<?(p), k(c) = fc 2 , /(c) = A 
and D c = 1, we obtain 

% = V-(Dh(p)Vp-xg(fi)Vc), (13) 



e—= Ac- k 2 c + Xp. (14) 

These equations are isomorphic to the nonlinear mean 
field Fokker-Planck equations ©-jS]) provided that we 
make the correspondance ^(r,i) = — c(r, t): the poten- 
tial of interaction is played by the concentration of the 
secreted chemical. 

It is important to note that the Kellcr-Scgel model is a 
mean field model. If we come back to the exact microscopic 
equations 

-±= X Vc d (r i (t),t) + V2DKi(t), (15) 



8 N 
e-^=Ac d -k 2 c d + \J2^-n(t)), (16) 

i=l 

governing the evolution of each particle (bacteria, cells, 
social insects,...), the mean field approximation leading to 
the standard Keller-Segel model l[T3] ) -([l"i| ) with h(p) = 1 
and g(p) = p amounts to neglecting statistical correla- 
tions; see |64I65I66I40I59I67I55| for details. Let us first 
assume that the system is Markovian and possesses no 
intrinsic memory in the sense that c<j(r,t) = — f u(\r — 

r'\)p d (r',t)dr' where p d (r,t) = £\ =1 6(r - r;(t)) is the 
exact distribution of particles. When the particles inter- 
act via a weak long-range binary potential u(|r — r'Q, it 

4 In this paper, we shall restrict ourselves to this situation. 
The original Keller-Segel model (TTT|) - ([12]) where Di and D 2 
depend on both p(r, t) and c(r, t) does not seem to possess 
the nice "thermodynamical" properties of the reduced Keller- 
Segel (|13[1 - (|14|1 such as Legendre structure of the free energy 
functionals, canonical //-theorem, Einstein relation etc. 
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can be shown that the mean field approximation becomes 
exact in a proper thermodynamic limit N — > +00 |40j . 
For example, in Eq. (|16[) these assumptions correspond to 
e = and k <C 1. By contrast, when memory is not zero 
and the interaction is short range, corresponding to e 7^ 
and fc 3> 1 in Eq. (|16[) . the mean field approximation can 
lead to wrong results. In the context of chemotaxis, the 
differences between mean field and non mean field models 
have been discussed by Grima [68] who showed the failure 
of the mean field approximation for the prediction of the 
width of the aggregate sizes (the disagreement becomes 
very severe close to the critical point where we know that 
mean field approximation breaks down in general) . On the 
other hand, the mean field approximation assumes that 
the number of particles N — > +00. In stellar systems and 
plasmas, for example, this is always the case. However, 
for biological systems, the number of interacting bacteria 
or cells is frequently less than a few thousands so it may 
be relevant to return to a microscopic description of the 
bacteria or cells' movement in terms of iV-body stochastic 
equations like (fT5| - (ri6)l as discussed in [67]. In this paper, 
however, we shall exclusively consider mean field models. 
As we have seen, mean field approximation works well for 
Markovian systems (e = 0) with weak long-range interac- 
tions if (i) N — > +00 and (ii) if we are not too close from a 
critical point. Mean field approximation breaks down: (i) 
for non-Markovian systems (e 7^ 0) (ii) close to a critical 
point (iii) for small values of N. 



2.5 Generalized free energy and H-theorem 

When <P ex t is an external potential, we define the energy 

by 



where C(p) is a convex function (C" > 0) defined by 



E = / p$ext dr. 



(17) 



When <!> is determined by Eq. the self- interaction en- 
ergy is given by 



E= - I P $dv. 



(18) 



Finally, when $ is determined by Eq. ([3]), the energy is 

E =^J [(V0) 2 + fc 2 <Z> 2 ] dr + J p$dr. (19) 

For e = 0, the expression (TT?j) reduces to Eq. (fT5|) . On the 
other hand, we define the temperature by 



D 

T =—. 

X 



(20) 



Therefore, the Einstein relation is preserved in the gener- 
alized thcrmodynamical framework. We also set f3 = 1/T. 
We introduce the generalized entropic functional 



S = - j C(p)dr, 



(21) 



C"{p) = 



Hp) 

9(PY 



(22) 



Note that the relation (f2"2"| defines the function C{p) up to 
a term of the form Ap + B. After integration over the do- 
main, the first term is proportional to the mass which is a 
conserved quantity and the second term is just a constant. 
Therefore, these terms play no physical role. However, we 
can adapt the values of the constants A and B in order to 
obtain convenient expressions of the entropy. 

Since the system is dissipative, the energy is not con- 
served. What is fixed instead of the energy is the temper- 
ature defined by the Einstein relation ([20]) . Since D cx T, 
the temperature measures the strength of the stochastic 
force in Eq. (p]). This corresponds to a canonical descrip- 
tion where the system is in contact with a heat bath. 
Note that the heat bath is completely disconnected from 
the long-range potential of interaction; it corresponds to 
short-range interactions modelled by the stochastic term 
in Eq. ([1]) . For an isolated system described by the micro- 
canonical ensemble the proper thcrmodynamical potential 
is the entropy. Alternatively, for a dissipative system de- 
scribed by the canonical ensemble, the relevant thermo- 
dynamical potential is the free energy. We introduce the 
generalized free energy 



F = E-TS. 



(23) 



The definition of the free energy (Legcndrc transform) is 
preserved in the generalized thcrmodynamical framework. 

When the energy is given by Eqs. (|17[) or (fT8|) . a straight- 
forward calculation (see Appendix [Gj shows that 



F 



1 



X9(p) 



(Dh(p)Vp + X9 (p)V<P) 2 dr. (24) 



When the energy is given by Eq. (fT9|) . we obtain 
F =-— I (A<P - k 2 <P - \pfdr 



1 



xg{p) 



{Dh(p)Wp + X g{p)V<P) 2 dr. 



(25) 



Therefore, F < 0. This forms an H theorem in the canoni- 
cal ensemble. It is sometimes useful to introduce the Massieu 
function 



J = S - pE, 



(26) 



which is related to the free energy by J — ~j3F. Clearly, 
we have J > 0. We can now consider particular cases: if 
D = (leading to T = 0), we get F = E so that E < 0. 
If x = (leading to = 0), we have J = S so that S > 0. 



2.6 Stationary solution 

The steady state of Eq. ((8]) satisfies F — leading to J = 
or explicitly 



Dh{p)v P + X g{p)v$ = 0. 



(27) 
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Using Eqs. (J20J) and ([22]). we get 

C"(p)Vp + /3V<2> = 0, 
which can be integrated into 

C'(p) = -0$-a, 



(28) 



(29) 



where a is a constant of integration. Since C is convex, 
this equation can be reversed to give 



p(r) = F(/?*(r)+a), 



(30) 



where F(x) = (C") _1 (— x) is a monotonically decreasing 
function. Thus, in the steady state, the density is a mono- 
tonically decreasing function p = p(<P) of the potential. 
We have the identity 



f/(0) = 



(3 



C"{p) 



(31) 



Substituting Eq. ([2]) in Eq. (|29|) . we find that the density 
profile is determined by an integro-diffcrential equation of 
the form 



C'{p) = -/? / p(r>(|r-r'|)dr' 



(32) 



More specifically, when the potential is given by Eq. ([JJ, 
we obtain a mean field equation of the form 



A<P-k 2 $ = \F(/3<P + a). 



(33) 



The constant of integration a is determined by the total 
mass M (which is a conserved quantity). Finally, we note 
that the generalized entropy |2T|) is related to the distri- 
bution bv [69]: 



C{p) 



F~ 1 (x)dx. 



(34) 



Equation (|30p determines the distribution from the en- 
tropy and Eq. (|34p determines the entropy from the dis- 
tribution. 



2.7 Minimum of free energy 

The critical points of free energy at fixed mass are deter- 
mined by the variational problem 



5F + Ta5M = 0, 



(35) 



where a is a Lagrange multiplier. We first consider the case 
where $ is given by Eq. ([2]). Therefore, the free energy that 
we consider is 



(36) 



F \P\ = \ J P^dv + T J C(p)dv. 



After straightforward calculations, we find that Eq. ([35 
leads to 



C'{p) = -[3$ -a. 



(37) 



Therefore, comparing with Eq. (j29p . we find that a sta- 
tionary solution of Eq. ([5]) is a critical point of F at fixed 
mass. Furthermore, it is shown in Appendix [B] that a 
steady state of Eq. ([8]) is linearly dynamically stable iff 
it is a minimum (at least local) of F at fixed mass. This 
property also results from Lyapunov's direct method [36] . 
Indeed, we have established that 



F < 0, 



F = *-> d t p = 0. 



(38) 



This implies that p(r) is linearly dynamically stable iff it 
is a (local) minimum of F at fixed mass (maxima or saddle 
points of F are dynamically unstable). In this sense, dy- 
namical and generalized thcrmodynamical stability in the 
canonical ensemble coincide. Furthermore, if F is bounded 
from below 0, we can conclude from Lyapunov's theory 
that the system will converge to a stable steady state for 
t — > +oo which is a (local) minimum of F[p] at fixed mass. 
If several local minima exist, the choice of the final steady 
state will depend on a complicated notion of basin of at- 
traction. 

In conclusion, a steady solution of the nonlinear mean 
field Fokker-Planck equation (J8]) is stable iff it satisfies the 
minimization problem: 



min 

p 



M[p] = M}. 



(39) 



Taking the second variations of F and using Eq. ([3Tj) . the 
condition of stability can be written 



>0, (40) 



for all perturbations 5p that conserve mass. If & ex t (r) is an 
external potential, the second term in Eq. (|40|) vanishes. 
Therefore, the second variations of the free energy are al- 
ways positive 5 2 F = -TS 2 S = (1/2)T/ C"(p){Sp) 2 > 
so that a critical point of F at fixed mass is necessarily a 
minimum. If the potential is given by Eq. Q , the free en- 
ergy is a functional F[p, <P] of the density p and potential 
(p. The cancellation of the first order variations of F with 
respect to 6p and yields Eqs. (|37[) and respectively. 
The condition of stability can be written 



5 2 F[6p,6$] 



{5pf 



dr + SpS<Pdr 



p>{$) 

--^J (AS® - k 2 S<P)S<Pdr > 0, (41) 

for all perturbations Sp and 8<P that conserve mass. From 
now on, we shall only consider the case where the potential 
is given by Eq. ([2]) since the case of Eq. ([3]) can be treated 
similarly by following the lines sketched above. 



5 There are important cases, like the system of self- 
gravitating Brownian particles, where the free energy is not 
bounded from below. In that case, the system can either relax 
towards a local minimum of F at fixed mass (when it exists) 
or collapse to a Dirac peak [56], leading to a divergence of the 
free energy F(t) — > — oo. 
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2.8 Particular cases 

If we take h(p) = 1 and g(p) = 1/C"(p), we get 



dp 
dt 



V • [ DVp- 



X 



C"(p) 



V<2> 



(42) 



In that case, we have a constant diffusion D{p) = D and a 
variable mobility \(p) = x/(pC"(p))- If we take g(p) = p 
and h(p) = pC"(p), we get 



dp 

at 



V • (DpC"(p)Vp + X pV$) 



(43) 



In that case, we have a constant mobility x{p) = X an d 
a variable diffusion .D(p) = Dp[C(p)/p]'. Note that the 
condition D(p) > requires that [C(p)/p]' > 0. This gives 
a constraint on the possible forms of C(p). 



2.9 Generalized Smoluchowski equation 

The NFP equation (|43|) can be put in the form of a gen- 
eralized Smoluchowski (GS) equation 



dp 
dt 



= V • [x(Vp + pV<£)] 



(44) 



with a barotropic equation of state p(p) given by 

p'{p) = T P C"{p). (45) 

Since C is convex, we have p'{p) > 0. On the other hand, 
integrating Eq. (|45|) twice, we get 



TC(p)=p[ Pp -^ d p'. 



(46) 



Therefore, the free energy ([56]) can be rewritten 

= ^ / P$dr + j pf yrdp'dr. (47) 

With these notations, the //-theorem (|24p becomes 
Y 



F= - I -(Vp + pV$) 2 dr < 0. 



(48) 



The stationary solutions of the GS equation ([4"4"|) satisfy 
the relation 



Vp + = 0, 



(49) 



which is similar to a condition of hydrostatic equilibrium. 
Since p = p(p), this relation can be integrated to give 
p = p(<P) through 



J p' 



This is equivalent to 



P'(p) _ 1 



p'm ' 



(50) 



(51) 



This relation can also be obtained from Eqs. (|4"5|) and 
(f3"Tj) . Therefore, we recover the fact that, in the steady 
state, p = p(^) is a monotonically decreasing function of 
<P. We also note the identity 



p{p) = \b{ P ) P = Tp 2 

A 



C( P ) 



= T[C'{p)p-C{p)}. 



(52) 



Finally, we note that the relation (|50|) . equivalent to 
the condition of hydrostatic equilibrium (|49|) . can also be 
obtained by extremizing the free energy (|47p at fixed mass 
writing 5F — a5M = 0. More precisely, we have the im- 
portant result: a steady solution of the generalized Smolu- 
chowski equation is linearly dynamically stable iff it 
is a (local) minimum of the free energy F[p] at fixed mass 
M[p] = M. This corresponds to the minimization problem 
p?]) with Eq. (gTJ). Note that the GS equation has been 
introduced here in a very general context. In this sense, 
it is valid for arbitrary value of x ( n °t necessarily small) . 
Alternatively, we shall see in Sec. 14.81 that a generalized 
Smoluchowski equation with \ = 1 — ^ can be derived 
from a generalized Kramers equation in the strong friction 
limit £ — > +oo. 



2.10 Phenomenological derivation of the nonlinear 
Fokker-Planck equation 

The form of the diffusion current appearing in the NFP 
equation (|8|) can be obtained by different phenomenolog- 
ical procedures. 



2.10.1 Functional derivative 

For a given free energy functional F[p], we can introduce 
phcnomcnologically a dynamical model by writing the evo- 
lution of the density as a continuity equation dtp = —V - J 
where the current is proportional to the gradient of the 
functional derivative of the free energy [36], i.e. 



dp 

dt 



= v 



SF 
Sp 



For the free energy ([36[) . we have 

SF 
dp 



so that 



^=V-fr(r,i) (TpC"(p)Vp + pV$)] 



(53) 



(54) 



(55) 



This equation is more general than Eq. ([S]). It shows that, 
for a given free energy, we can introduce an infinite class 
of NFP equations where x(r, t) is an arbitrary positive 
function of position and time. In particular, %(r, t) can be 
a positive function of the density p(r, i). If we set x(r, t) = 



P.H. Chavanis: Nonlinear mean field Fokker-Planck equations. 



9 



Xd(p)/p we recover Eq. ([8]) as a particular case. We can 
also write Eq. (|55p in the alternative form 



dp 
di 



X(r,i) \TVp- 



1 



-v<z> 



(56) 



C"(p) 

where x(r, t) is an arbitrary positive function of position 
and time. If we set x(r,£) = xM/ 9 ) we recover Eq. (J8|). 
These two alternative forms (|55l) and (|56|) were given in 
03]. 

On the other hand, the formulation (|53|) ensures that 
the free energy decreases monotonically provided that x(r, t) 
is positive. Indeed, 



F 



6p 



/• 



5Fdp , 

- P m dv = - 

X (r, V — 



SF 

2 



dr < 0. (57) 



Furthermore, a steady state of Eq. (|53|) satisfies f = 0, i.e 
S7(SF/5p) = leading to 

6 -f = -To, (58) 

where a is a constant of integration. This is equivalent to 
the extremization of F at fixed mass M if we write the 
first variations as 5F + Ta6M = 0. Therefore, a steady 
state extrcmizes the free energy at fixed mass. Using Eq. 
(|54p , we find that the equilibrium density profile satisfies 



C'(p) = -(3$ -a. 



(59) 



Finally, using Lyapunov's direct method, one can show 
that a steady state p(r) of Eq. ([53")) is linearly dynami- 
cally stable iff it is a (local) minimum of F at fixed mass. 
We again emphasize that the stationary states of Eq. ([55]) , 
and the iJ-theorem, only depend on the form of the free 
energy F[p] and are independent on the positive function 
x(r, t). These properties are therefore valid for the whole 
class of NFP equations associated with a given free en- 
ergy functional. Therefore, for a given free energy specified 
by C(p), we can construct an infinite class of NFP equa- 
tions with arbitrary x(r, t) that have the same equilibrium 
states (|59|) but a different dynamics. 

2.10.2 Onsager's linear thermodynamics 

The previous approach is equivalent to Onsager's linear 
thermodynamics. Indeed, at equilibrium, we expect that 
the distribution p(r) minimizes the free energy F at fixed 
mass. This leads to Eq. (|58|) or (|59|) . Noting that the chem- 
ical potential 

\{T,t) = ?f=TC(p)+$, (60) 
dp 

is uniform at equilibrium, the linear thermodynamics of 
Onsager suggests that, close to equilibrium, the current is 
proportional to the gradient of the chemical potential. If 
we write 

J = -xM)pVA(r,i), (61) 
the linear thermodynamics of Onsager leads to Eq. (|53| . 



2.10.3 Maximum Free Energy Dissipation Principle 

The same results can be obtained from a variational for- 
mulation, called the Maximum Free Energy Dissipation 
(MFED) principle which is the canonical ensemble ver- 
sion of the Maximum Entropy Production (MEP) prin- 
ciple [33] . At equilibrium, the optimal distribution p(r) 
minimizes the free energy F[p] at fixed mass M. Out of 
equilibrium, we may expect that the optimal current J 
maximizes the rate of free energy dissipation F[3] under 
some constraints. This can be viewed as a variational for- 
mulation of Onsager's linear thermodynamics. The rate of 
dissipation of free energy is given by 



F = 



6F dp 
Jp~~dt 



dr=- 



5F_ 
Jp~ 



V • J dr = 



dp 



(62) 



We shall determine the optimal current J* which maxi- 
mizes the rate of dissipation of free energy F under the 
constraint J 2 < C(r, t) putting a physical bound on |J|. It 
can be shown that the bound is always reached so that we 
can replace the inequality by an equality. Thus, we write 
the variational problem as 



SF + 5 



dr ) = 0, 



(63) 



where x( r )^) is a local Lagrange multiplier. Performing 
the variations on J, we obtain 



J* = -X {r,t)pV— , 
dp 



(64) 



which returns Eq. ([53)1 . Note that if we introduce the "dis- 
sipation" function 



E,, 



2px(r,t) 



■ dr, 



we have 



(65) 



(66) 



On the other hand, since 5 2 {F + E d ) = - J ^~dr < 0, 



the optimal current (|64p maximizes the dissipation of free 
energy under the constraint J 2 < C(r, t). 



2.11 Kinetic derivation of the nonlinear Fokker-Planck 
equation 

Nonlinear Fokker-Planck equation where the diffusion co- 
efficient and the mobility explicitly depend on the local 
concentration of particles can be derived from a kinetic 
theory, starting from the master equation, and assuming 
that the probabilities of transition explicitly depend on the 
occupation numbers (concentrations) of the initial and ar- 
rival states. We briefly summarize the approach developed 
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by Kaniadakis |12| and make the link with the phenomeno- 
logical equations studied previously. 

We introduce a stochastic dynamics by defining the 
probability of transition of a particle from position r to 
position r'. Following Kaniadakis [Hj, we assume the fol- 
lowing form 



7r(r 



w(r,r-r')a[p(r,t)]b[p(r',t)}. (67) 



Usual stochastic processes correspond to a(p) = p and 
b{p) = 1: the probability of transition is proportional to 
the density of the initial state and independent on the den- 
sity of the final state. They lead to the ordinary Fokker- 
Planck equation ([86]) as will be shown below. Here, we 
assume a more general dependence on the occupancy in 
the initial and arrival states. This can account for micro- 
scopic constraints like close-packing effects that can in- 
hibitate the transition. Quite generally, the evolution of 
the density satisfies the master equation 



dp 
dt 



[ir(r' -*•!•)- 7r(r -> r')]dr'. 



(68) 



Assuming that the evolution is sufficiently slow, and local, 
such that the dynamics only permits values of r' close to r, 
one can develop the term in brackets in Eq. (|68[) in powers 
of r — r'. Proceeding along the lines of [T2], we obtain a 
Fokkcr-Planck-like equation 



dp d 

dt dxi 

with 



c* + §f ) 



din k(p) dp 
dp %3 dxj 



lip) = a{p)b(p), n{p) = 



ajp) 



and 



yiw(r,y)dy, 



(69) 
(70) 

(71) 



Cij'W = I; J yiyjw(r,y)dy. (72) 

The moments Q and Qj are fixed by the ordinary Langevin 
equation 



dr 



- = - X V$ + V2DTL(t), 



(73) 



where % and D are constant. Assuming isotropy Q = Ji 
(ij = DSij, the kinetic equation (|69|) becomes 



dp 
dt 



V 



(J + V^) 7 (p)+ 7 (p)^^^Vp 
dp 



(74) 



Now, according to the Langevin equation (|73p . D is inde- 
pendent on r and J = %V<?. Thus, we get 



dt 



^7(P)^^V P + X7(P)V<P 



If we define 

Hp) = lip) 



dp 

d\nn(p) 
dp 



g{p) = lip), 



(75) 



(76) 



the foregoing equation can be written 
dp 



dt 



v ■ [Dh{fi)V P + xg(p)v$] 



(77) 



and it coincides with the phenomenological equation JS}. 
We note that 

\n K (p) = C'(p). (78) 
We also have the relations 

<P) = VliP)<P) = VW)e C ' ip)/ '\ (79) 



Kp) 



1 jip) 

k( P ) 



Vaip) 



-C"(p)/2 



Inversely 

g(p) = a(p)b(p), C'ip) = In 

h(p)=b(p)a'(p)~a(p)b'(p). 



aip) 



Hp) 



(80) 



(81) 



(82) 

It seems natural to assume that the transition proba- 
bility is proportional to the density of the initial state so 
that a{p) = p. In that case, we obtain an equation of the 
form 



dp 
dt 



V • (D [b(p) - pb'(p)) Vp + xpbipm) 



(83) 



Note that the coefficients of diffusion and mobility are not 
independent since they are both expressed in terms of 6(p). 
Choosing b[p) = 1, i.e. a probability of transition which 
docs not depend on the population of the arrival state, 
leads to the standard Fokker-Planck equation (f86|) . If, 
now, we assume that the transition probability is blocked 
(inhibited) if the concentration of the arrival state is equal 
to (To, then it seems natural to take b{p) = 1 — p/<tq. In 
that case, we obtain 



^ = V • (£Vp + X p(l - PM))V<2>) 



(84) 



which will be considered in Sec. 13.51 Inversely, we can 
wonder what the general form of the mobility will be if 
we assume a normal diffusion h(p) = 1. This leads to 
b(p) — pb'ip) — 1 which is integrated in b(p) = 1 + Kp 
where if is a constant. Interestingly we find that this con- 
dition selects the class of fermions (K = — 1) and bosons 
(K = +1) and intermediate statistics (arbitrary K). The 
corresponding NFP equation is 



dp 
dt 



V ■ (DVp + xpil + Kp)V$) 



which will be considered in Sec. 



(85) 



In Sec. 12.11 we have obtained generalized Fokker-Planck 
equations by using ordinary Master equations (based on usual 
transition probabilities a(p) — p and b(p) = 1) and general- 
ized Langevin equations where the diffusion coefficient and the 
mobility depend on the density. In this section, we have ob- 
tained generalized Fokker-Planck equations by using general- 
ized Master equations (based on density dependent transition 
probabilities) and ordinary Langevin equations with constant 
coefficients. 



P.H. Chavanis: Nonlinear mean field Fokker-Planck equations. 



11 



3 Examples of nonlinear Smoluchowski 
equations and generalized Keller-Segel 
models 

In this section, we give several explicit examples of nonlin- 
ear mean field Fokker-Planck equations. Some correspond 
to well-known forms of entropies, and others are new. We 
emphasize that these equations can have applications in 
different domains of physics as discussed in Sec. 12.31 Im- 
portantly, they can provide generalized Keller-Segel mod- 
els of chemotaxis (see Sec. l2.4[) . Most of these models have 
not been considered before in biology because the connec- 
tion with generalized thermodynamics was not made. This 
is why we give a relatively detailed description of these 
models since their applications in biology are new. 



3.1 Standard model: Boltzmann entropy 

If we take h(p) — 1 and g(p) = p, we get the ordinary 
Smoluchowski equation 



|£ = V.(I>V/> + xpV#). 



(86) 



It corresponds to an ordinary diffusion D(p) = D and 
a constant mobility x(p) = X- The associated stochastic 
process is 



dr 
~dt 



The entropy is the Boltzmann entropy 



S 



plnpdr, 



and the stationary solution of Eq. 
distribution 



(87) 



(88) 



is the Boltzmann 



-fi$—a—l 



The pressure law is 



p(p) = P T - 



(89) 



(90) 



This is similar to the equation of state for an isother- 
mal gas with constant temperature T . When the Fokker- 
Planck equation l[8"6"|) is coupled to the Poisson equation 
([5]), we obtain the Smoluchowski-Poisson system describ- 
ing a gas of self-gravitating Brownian particles 
When the Fokker-Planck equation ([86]) is coupled to the 
field Eq. ([3]), we obtain the standard Keller-Segel model 
describing the chemotactic aggregation of biological pop- 
ulations l52l. 



It corresponds to a power law diffusion D(p) = Dp 11 " 1 
and a constant mobility x(p) — X- This equation was 
introduced by Plastino & Plastino [21]. The associated 
stochastic process, introduced by Borland [25], is 



dr 

It 



— = -xV0 



-n(t). 



(92) 



This model can take into account effects of non-ergodicity 
and nonextensivity. It leads to a situation of anomalous 
diffusion related to the Tsallis statistics. For q = 1, we 
recover the standard Brownian model with a constant dif- 
fusion coefficient, corresponding to a pure random walk. 
In that case, the sizes of the random kicks are uniform 
and do not depend on where the particle happens to be. 
For q ^ 1, the size of the random kicks changes, depend- 
ing on the distribution of the particles around the "test" 
particle. A particle which is in a region that is highly pop- 
ulated [large p(r, t)] will tend to have larger kicks if q > 1 
and smaller kicks if q < 1 . Since the microscopies depends 
on the actual density in space, this creates a bias in the 
ergodic behavior of the system. Then, the dynamics has 
a fractal or multi- fractal phase space structure [25]. The 
generalized entropy associated to Eq. ([9Tjl is the Tsallis 
entropy 

8= l —({ P «-p)dv 1 (93) 

3- 1 J 

and the stationary solution is the Tsallis distribution 

i 

'~ 1 [l-(«-l)03# + a)£ /(, - 1) . (94) 



The pressure law is 



p(p) = Tp q . 



(95) 



This is similar to a polytropic gas with an equation of 
state p = Kp 1 where K = T plays the role of a polytropic 
temperature and q = 7 is the polytropic index (we also 
set 7 = 1 + 1/n). Note that the Tsallis entropy can be 
written 



8 = - J P^(g) Pdr, 

re we have introduced the g-logarithm 
1 



ln (g) (x) 



-(x^ 1 - 1). 



9-1 

The stationary solution can be written 



(96) 



(97) 



3.2 Power law diffusion: Tsallis entropy 

If we take h(p) = qp q ~ x and g(p) = p, we obtain the NFP 
equation 



dp 



V • (DVp q + x/?V<£) 



(91) 



P = 



1 \ "~ 



(?) 



with the g-exponential 

e {q) (x) = [l + {q-l)x]^. 



(98) 



(99) 
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For q = 1, we recover the standard model (|86|) . For q = 
2, we have some simplifications. In that case, the NFP 
equation (f9Tj) becomes 



dp 

at 



= V • (D\7p 2 + x/?V0) 



The entropy is the quadratic functional 



P dr, 



and the stationary solution is 

p= -l(/3$ + a), 



(100) 



(101) 



(102) 



corresponding to a linear relation between the density and 
the potential. In that case, the differential equation (I33|) 
determining the steady state reduces to the Helmholtz 
equation. Finally, the pressure is 



p(p) = Tp 2 



(103) 



corresponding to a polytrope with index n — 1. When 
the NFP equation (|9Tj) is coupled to the Poisson equation 
(fSJ), we obtain the polytropic Smoluchowski-Poisson sys- 
tem describing self-gravitating Langevin particles. When 
the NFP equation dUD is coupled to the field Eq. ©, we 
obtain a generalized Kellcr-Scgcl model of chemotaxis tak- 
ing into account anomalous diffusion. These models have 
been introduced by Chavanis & Sire [45151153] , 



3.4 Power law diffusion and drift: Tsallis entropy 

We introduce here a new model generalizing the polytropic 
model ([£[]). If we take h(p) = qp^^- 1 and g(p) = 
we obtain 



dp 
lit 



(108) 



This corresponds to a power law diffusion D(p) 
and a power law mobility \(p) = \p^ . The associated 
stochastic process is 



dr 

~dt 




2D ±p^R(t). (109) 



Since p^ can be put in factor of the diffusion current in 
Eq. (|108p . this model belongs to the infinite family of NFP 
equations associated to the Tsallis entropy with index q 
(see discussion in Sec. 12. 10~T)l . 

For p = 0, we recover Eq. (|91[) with a constant mobility 
and a power law diffusion. For (p,q) = (0,0), we recover 
the logotropic Smoluchowski equation (|104|) provided that 
we make the transformation Dq — > D. For p = 1 — q, we 
have a normal diffusion and a power law mobility 



^ = V.(DqVp + X p 2 - q V<P) 



(110) 



3.3 Logotropic distributions: log-entropy 

If we take h(p) = 1/p and g(p) = p, we obtain 

dp 



at 



V ■ (DV In p + xpV0) 



(104) 



The generalized entropy associated to Eq. (|104[) is the log- 
entropy 



lnpdr, 



and the stationary solution is 



1 



(105) 



(106) 



For a quadratic potential <P e 



2 /2, this corresponds to 



the Lorcntzian function. The pressure is 

p(p) = T\np. (107) 

This is similar to a logotropic equation of state [7D] . This is 
also connected to a polytropic equation of state (or Tsallis 
distribution) with 7 = q = 0. Indeed, the logotropic model 
(|104|) can be deduced from Eq. (|9"Tj) by writing D\7p q = 
Dqp q ^ 1 W p, taking q = and rc-defining Dq — > D. When 
the NFP equation (|104|) is coupled to the Poisson equa- 
tion ([5]), we obtain the logotropic Smoluchowski-Poisson 
system. When the NFP equation (I104|) is coupled to the 
field Eq. ([3]), we obtain a generalized Kellcr-Scgel model of 
chemotaxis. These models have been introduced by Cha- 
vanis & Sire [71] . 



For q = 2, we get 

dp 



^=V-(2£V P + X V<?>), 



(Ill) 



which has the same equilibrium states as Eq. (|100[) . If 
we assume furthermore that is given by the Poisson 
equation ([5]), Eq. (|111[) reduces to the linear equation 



dp 

at 



IDAp + x Xp- 



Finally, for q = (making the transformation qD 
we obtain 

^ = V-0DVp + X p 2 Vf), 



(112) 



(113) 



which has the same equilibrium states as Eq. (|104| . When 
the NFP equation (|108|) is coupled to the field equation 
([3]), we obtain a generalized Kcller-Segel model of chemo- 
taxis taking into account anomalous diffusion and anoma- 
lous mobility. This model will be studied in a forthcoming 
paper, in continuity with [45j . 



3.5 Filling factor: Fermi-Dirac entropy 

If we take h(p) = 1 and g(p) = p(l — p/ctq), wc obtain 
dp 



dt 



V • (DVp + XP0- - p/^o)V^) 



(114) 
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This corresponds to a normal diffusion D(p) = D and 
a mobility x(p) = xO- ~ p/ a o) vanishing linearly when 
the density reaches the maximum value p max = ao- The 
associated stochastic process is 



dr 



■x(l-p/ao)V$ + V2DTl(t). 



(115) 



The generalized entropy associated with Eq. (|114[) is the 
Fermi-Dirac entropy in position space 



S 



-00 



— In — 

e>o a 



1 - 



CO 



In 1 - 



CO 



dr, 
(116) 



positive function x( r i^) m order to obtain a more general 
model with the same entropy and the same equilibrium 
states. Note finally that an equation similar to Eq. (|114[) 
has been introduced by Robert & Sommeria [ST] (see also 
28J) in the statistical mechanics of two-dimensional tur- 
bulence for two vorticity levels and ao- In that case, p 
represents the coarse-grained vorticity ZU and <P plays the 
role of the stream function ip. The "exclusion principle" 
leading to the Fermi-Dirac entropy (|116[) is a consequence 
of the 2D Eulcr equation implying that the vorticity levels 
cannot overlap so that w(r,t) < ao- These analogies be- 
tween chemotaxis and 2D turbulence are further discussed 
in [731. 



and the stationary solution is the Fermi-Dirac distribution 
in position space 



P 



<70 



1 + e ^+ a ' 



(117) 



From Eq. (jTTTJ) , we see that, in the stationary state, p < 
<tq. This bound is similar to the Pauli exclusion principle 
in quantum mechanics. In fact, we can show that p(r, t) 
remains bounded by <to during the whole evolution. For 
<To - ► +oo, we recover the standard model (|86|) . 

An alternative model, with the same entropy and the 
same equilibrium states, is obtained by taking h{p) = 
1/(1 — p/ctq) and g(p) = p. This leads to 



dp 
dt 



V • (-Da V ln(l - p/a ) + XP V$) . (118) 



This corresponds to a nonlinear diffusion with D(p) = 
—ao(D/p) ln(l — p/(T ) and a constant mobility \{p) = X- 
Equation (|118[) can be put in the form of a generalized 
Smoluchowski equation (|44[) with a pressure law 



(119) 



p(p) = -To- O ln(l - p/a Q ). 



For p <C ao, we recover the "isothermal" equation of state 
p = pT leading to the standard model (|86|) . However, for 
higher densities, the equation of state is modified and the 
pressure diverges when p — > (To ■ This prevents the density 
from exceeding the maximum value do- 
In the context of chemotaxis, the model (|114|) has 
been introduced by Hillen & Painter [72] and, indepen- 
dently, by Chavanis [33166173] , It provides a regulariza- 
tion of the standard Kcller-Segel model preventing over- 
crowding, blow-up and unphysical singularities. The filling 
factor (1 — p/<tq) takes into account the fact that the par- 
ticles cannot interpenetrate because of their finite size a. 
Therefore, the maximum allowable density is ao ~ l/a d . 
It is achieved when all the cells are packed together. In the 
model (|114p , it is assumed that the mobility vanishes when 
the density reaches the close packing value (p — > (To) while 
the diffusion is not affected. The alternative model ([118D 
has been introduced in Chavanis |33l73j . In that case, the 
mobility is assumed to be constant and the regularization 
preventing overcrowding is taken into account in the pres- 
sure law (|119|) . As explained in Sec. 12.10.11 we can multi- 
ply the diffusion term and the mobility term by the same 



3.6 Fermi, Bose and intermediate statistics 

If we take h(p) = 1 and g{p) = p(l + Kp), we obtain 

% = V • {DVp + xp(l + Kp)V$) . (120) 
at 

This corresponds to a normal diffusion D{p) = D and 
a variable mobility x{p) = x(l + Kp). The associated 
stochastic process is 



dr 

dt 



— = -x(l + Kp)V$ + V2DK(t). 



(121) 



Here, K is a real number taking positive or negative val- 
ues. When K > the mobility is enhanced in regions of 
large densities and when K < 0, it is reduced. This takes 
into account inclusion (K > 0) or exclusion (K < 0) prin- 
ciples. For K = 0, we recover the standard model (j86|) . 
The generalized entropy associated with Eq. (j!20p is 



plnp 



Kp) ]n(l + Kp) 



dr, 



and the stationary solution is 



1 



J4>+a _ K 



(122) 



(123) 



For K = +1 we obtain the Bosc-Einstcin statistics and for 
K = —1 we obtain the Fermi-Dirac statistics. For other 
values of K , we obtain intermediate statistics (quons) in- 
terpolating between fermions and bosons. For K = 0, we 
recover the Boltzmann statistics. 

An alternative model with the same entropy and the 
same equilibrium states is obtained by taking h(p) = 1/(1+ 
Kp) and g(p) = p. This leads to 



dt 



lvin(l 



Kp) + xpV<£ 



(124) 



This corresponds to a nonlinear diffusion such that D(p) = 
(D/Kp) \n(l + Kp) and a constant mobility x(p) = X- The 
pressure law is 



p(p) = — +Kp). 



(125) 
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We recall that Eq. (|120|) can be obtained from the master 
equation (|68|) when the transition probabilities are of the 
form (f67|) with a(p) = p and b(p) = 1 + Kp. Alternatively, 
the model ()124j) corresponds to a(p) 

Kp) 



pj\J\ + Kp and 



y/1 + Kp. 

The NFP equations <pT20|) and (fT24|) have been intro- 
duced by Kaniadakis & Quarati [27] (see also [28] in the 
context of the violent relaxation of 2D vortices and stel- 
lar systems). For K = +1, they can provide a dynamical 
model of the Bose-Einstein condensation in the canonical 
ensemble which has been studied in detail in [30] . When 
coupled to the field equation ^ , the NFP equations (|120p 
and (|124|) could also provide generalized Kcllcr-Scgel mod- 
els of chemotaxis. 



3.7 Mixed model: anomalous diffusion and filling factor 

The previous models focus individually on two important 
effects: anomalous diffusion (see Sees. I3.2ll3.4j) and exclu- 
sion constraints when the density becomes too large (see 
Sec. I3.5p . Here we introduce a mixed model which com- 
bines these two effects in a single equation. If we take 
h{p) = qp q+ ^ x and g(p) = p' I+1 (l — p/ct ), we obtain 

|£ = V • (Dqpfi+P-^p + XP» +1 (1 - p/ct )V#) . 



(126) 

This corresponds to a power law diffusion such that D(p) = 
[Dq/{q + m)]p 9+a '~ 1 and a mobility X (p) = XP M (1 - p/ao). 
The associated stochastic process is 



dr ... / 2Dg g+n-i _ , , 

-n = -XP fl (l - p/°o)V<P + J — ip->-R t). 
at v q + p 



(127) 



The generalized entropy corresponding to Eq. (|126|) is ob- 
tained by integrating twice the relation 



C"(p) 



A first integration gives 



qp 



,9-2 



1 - p/ao 



C'(p) = qa^-i U-\ , (129) 



where 



-dx. 



I 1-x 

Therefore, the generalized entropy can be expressed as 



(130) 



p/cro 

C(p)=qo q l *,_2 (*)<**• 



(131) 



Note that it does not depend on p since the term p M can 
be put in factor of the diffusion current in Eq. (|126p : see 
the discussion in Sec. 12.10.11 



Let us consider some particular cases, (i) For q = 1, 
Eq. (|126p has the same entropy and the same equilibrium 
states as Eq. (|114|) . (ii) For cto — > +oo, we recover Eq. 
(|108j) . (iii) For p = and g = 2, we have 



(132) 



^ = V • (£Vp 2 + X p(l - p/ct )V<P) . 
The generalized entropy is 



S = -2ai 



CTQ 



l_ZWfl-^-|dr. 



CT 



and the stationary solution is 



cto 



_ (/3<f+Q)/2<T 



(133) 



(134) 



For cto — > +oo, we recover Eq. (|102[) . We can also consider 
the alternative model 



dp 

dt 



2pD 



1 - p/ctq 



Vp + xpV<? , 



(135) 



which has the same entropy and the same equilibrium 
states as Eq. (| 132|) . The pressure law is 



P(p) = - 2T(J o M 1 - P/ a o) - p/ctq] • 



(136) 



(iv) For (p, q) — (0, 0) and performing the transformation 
qD — ► D, or directly taking h(p) = 1/p and g{p) = p(l — 
P/cto), we obtain 

^ = V- (.DVlnp + xp(l-p/CTo)V<£). (137) 

This corresponds to a logarithmic diffusion and a mod- 
ified mobility taking into account an exclusion principle 
through the filling factor. The generalized entropy is ob- 
tained from the relation 



C"(p) 



1 



P 2 (1-p/cto)' 



(128) leading to 



C' ( p)=-J-{ln(^ 

CT0 I VP 



and finally to the explicit expression 



1 + — 



S 



1 



In 



CTO 



cto/ VP 
We can consider the alternative model 

D 



1 dr. 



dt 



p(l - p/ct ) 



Vp + xpv<? 



(138) 



(139) 



(140) 



(141) 



with the same entropy and the same equilibrium states. 
The associated pressure law is 



P(P) 



-Tin 



ctq 



1 



(142) 
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3.8 K-entropy 

Wc consider the K-entropy 



(p 1+K - p 1 - K )dr. 



(143) 



This entropy was introduced by Kaniadakis [12| . It can be 
written 



(144) 



S K = - J /oln (K) pdr, 
with the K-logarithm 

= l-(x H -x-"). 



We have 



(145) 



C"(P) = ^[(1 + «)P K -(1- K )P"1, (146) 



with a = k + r and b = n — r. Of course, we could give 
many other examples of generalized Fokker-Planck equa- 
tions since there exists an infinite number of distributions 
and entropic functionals. Therefore, we found it more con- 
venient in [33] to formulate the problem in a general set- 
ting, using an arbitrary entropic functional of the form 

(ED. 



4 Nonlinear mean field Fokker-Planck 
equations in phase space 

We now describe nonlinear mean field Fokker-Planck equa- 
tions in phase space taking into account the inertia of the 
particles. Overdamped models will be recovered in a limit 
of strong friction. 



4.1 Generalized Kramers equation 



C"(p),-[(l + K )p K +(l-^1. (147) 

If we take g(p) = p and h(p) = pC"(p), we obtain the 
NFP equation 



dp _ /D 

i = v -U v(p 



+ K _i_ 1 — K 



) + XPV<£ 



(148) 



This corresponds to a power law diffusion D(p) = + 
p~ K ) and a constant mobility x{p) = X- The associated 
stochastic process is 



dr 

dt 



^-xV'P + VD^ + p-^/z-Rit), 



(149) 



and the stationary solution of Eq. (|148[) can be written 

i 

1 — K 



P = 



1 + K 

with the ^-exponential 



-(/3*+Q)/ v / r r ^ 5 



e (K) (x) = (kx + \/T+ k 2 x 2 ) 1/k . 
Finally, the pressure law is 

p(p) = \t(p 1+k + p 1 -"). 

For k = 0, we recover the standard model 
results can be generalized to the («, r) entropy 



Sk.t — 



(150) 

(151) 

(152) 
These 

(153) 



which reduces in some special cases to the Tsallis [9], 
Abe [TU] and Kaniadakis [T2] entropies. The corresponding 
NFP equation can be written 



dt 



l+a 



bp 



l-b\ 



, (154) 



Wc consider a system of N particles in interaction whose 
dynamics is described by the stochastic Ito-Langevin equa- 
tions 



dt 



Vi, 



(155) 



dvi 
~dt 



-f (/i)vi - V# 4 + ^2D(fi)TLi(t), (156) 



where <£(r, t) is a self-consistent potential given by the 
mean field Eq. @. In ordinary models, the friction £ and 
the diffusion coefficient D are constant. In that case, the 
statistical equilibrium state is the Maxwcll-Boltzmann dis- 
tribution / ~ e~° e where e = v 2 /2 + ^(r) is the individ- 
ual energy and the temperature T = 1/(3 is given by the 
Einstein relation T = D/t;. Here, for sake of generality, 
the friction coefficient £(/) and the diffusion coefficient 
D(f) are allowed to depend on the distribution function 

/(r,v,i) = (Y^iS(r -T t (t))5{v -Vi{t))). This can take 
into account microscopic constraints that affect the dy- 
namics and modify the equilibrium distribution. The evo- 
lution of the distribution function /(r, v, t) is governed by 
the nonlinear mean field Fokker-Planck equation 



df df „ x df d 

— + v • — - W ■ — = — 
dt dr dv dv 



d_ 



(D(f)f) +f(/)/v 



(157) 



coupled to Eq. $2§. We introduce the notations 



Dh(f) = -n{fD{f)), 
df 



&(/) = /£(/), (158) 



where D and £ are positive constants and h(f) and g(f) 
are positive functions. The ordinary model with constant 
diffusion D(f) = D and constant friction £(/) = £ is re- 
covered for h(f) = 1 and g(f) = f. With these notations, 
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the NFP equation (|157|) can be rewritten in the form of a 
generalized Kramers (GK) equation 



It can be put in the conservative form 



where 



df_ _9_ 

dt dv ' 



Dh(f)^+^g(f)v 



(159) 
(160) 
(161) 



is a diffusion current. This structure guarantees the con- 
servation of mass M = J fdrdv. 



4.2 Generalized free energy and H-theorem 

We define the energy by 

E =\J A' 2 drdv + 5 J P®dr = K + W, (162) 

where K is the kinetic energy and W is the potential en- 
ergy. We define the temperature by 



T 



D 

T 



(163) 



The Einstein relation is preserved in the generalized ther- 
modynamical framework. We introduce the generalized 
entropic functional 



C(f)drdv, 

where C(f) is a convex function defined by 

Hf) 



c"(f) 



9(f)' 



(164) 



(165) 



Finally, we introduce the generalized free energy 

F = E-TS. (166) 

The definition of the free energy (Lcgcndre transform) is 
preserved in the generalized thermodynamical framework. 
Explicitly, 

F[f] = \J fv 2 drdv + lJ p <Pdr + Tj C(f) drdv. 

(167) 

A straightforward calculation (see Appendix [G| shows 
that 

2 



Therefore, F < 0. We can also introduce the Massieu func- 
tion J = S — f3E which is related to the free energy by 
J = -0F. It satisfies j > 0. If D = (leading to T = 0), 
we get F = E so that E < 0. If % = (leading to = 0), 
we get J = S so that S > 0. 



4.3 Stationary solution 



The steady states of Eq. (|159[) must satisfy F = 0. Ac- 
cording to Eq. (|168j) , this implies J = or explicitly 



Dh(f)?f+Zg(f)v 
ov 



Using Eqs. (fT63|) and ifTBg]) . we get 



C"(f) 



dv 



0v = 0, 



which can be integrated into 



C'(f) = -0 



A(r) 



(169) 



(170) 



(171) 



where A(r) is a function of the position. Since C is convex, 
this relation can be reversed to give 



/(r,v)=F 



/?( y + A(r 



(172) 



where F{x) = (C')~ 1 (—x) is a decreasing function. Since 
J = and df /dt = 0, the steady solution of Eq. (|159[) 
must also cancel the advective term 



df „ x df 
-J- _ ■ — 

dv dv 



0. 



(173) 



In other words, the steady solution of Eq. (|159|) is a par- 
ticular stationary solution of the Vlasov equation (l.h.s.) 
whose form is selected by the "collision" term (r.h.s.). Sub- 
stituting Eq. (fT72"l) in Eq. ([173]) . we get 



(VA- V<P) - v = 0, 
which must be true for all v. This yields 
A(r) = <P(r) + a/ (3, 



(174) 



(175) 



where a is a constant of integration. Therefore, the sta- 
tionary solution of Eq. (| 1 59[) is given by 



C'{f) = -f3e 



or, equivalcntly, 



/(r,v) = F[(3e(r,v) + a], 



(176) 



(177) 



where e = v 2 /2 + <P(r) is the energy of a particle. Thus, 
in the steady state, the distribution function / = /(e) 
is a function of the individual energy The potential <& is 
determined by an integrodifferential equation obtained by 
substituting Eq. (flTT]) in Eq. using p = J fdv. The 
constant a is determined by the conservation of mass. On 
the other hand, differentiating Eq. (|176|) . we obtain 



df_ 



c"(jy 



(178) 



Since C is convex, i.e. C" > 0, the preceding relation 
implies that /'(e) < 0. Therefore, /(e) is a decreasing 
function of the energy. 
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4.4 Minimum of free energy 



4.6 Functional derivative 



The critical points of free energy at fixed mass are deter- 
mined by the variational problem 



5F + Ta5M = 0, 



(179) 



where a is a Lagrange multiplier. These variations give 



C'(f) = -f3e - a. 



(180) 



Therefore, comparing with Eq. (|176| . we find that a sta- 
tionary solution of Eq. (|159|) is a critical point of F at 
fixed mass. Furthermore, it is shown in [33] that a steady 
state of Eq. (|159j) is linearly dynamically stable iff it is a 
minimum (at least local) of F at fixed mass. In this sense, 
dynamical and generalized thermodynamical stability in 
the canonical ensemble coincide. This property also re- 
sults from Lyapunov's direct method [33]. Finally, if F is 
bounded from below, we conclude from the above proper- 
ties that the system will converge to a stable steady state 
for t — > +oo which is a (local) minimum of F[f] at fixed 
mass. If several local minima exist, the choice of the fi- 
nal steady state will depend on a complicated notion of 
basin of attraction. In conclusion, we have the important 
result: a steady solution of the generalized Kramers equa- 
tion \159\) is linearly dynamically stable iff it is a (local) 
minimum of the free energy F[f] at fixed mass M[f] = hi . 
This corresponds to the minimization problem: 



nun {F[f} | M[f]=M}. 



(181) 



Taking the second variations of F and using Eq. Q178[) . 
the condition of dynamical stability is 

S 2 F[Sf] =-\[J jZ drdv - J S P S<Pdr } * °> ( 182 ) 
for all perturbations 8f that conserve mass. 

4.5 Particular cases 

If we take h(f) = 1 and g(f) = 1/C"{f) we get 

df df „ x df d ( df £ 

— + v • — - V<£ ■ — = — • D— A — v 

dt dv dv dv \ dv C"(f) 

(183) 

In that case, we have a constant diffusion D(f) = D and 
a variable friction £(/) = £/[/C"(/)]. If we take g(f) = f 
and h(f) = fC"(f) we get 

at or ov ov \ ov 

(184) 

In that case, we have a constant friction £(/) = £ and a 
variable diffusion D(f) = Df[C(f)/f}'. 



For a given free energy functional F[f], we can introduce 
phcnomcnologically a dynamical model by writing the evo- 
lution of the distribution function as a continuity equation 
dtf = — • J where the current is proportional to the gra- 
dient in velocity space of the functional derivative of the 
free energy, i.e. 



df df df d 

— + v — - W ■ — = — 
dt dr dv dv 



For the free energy (|167[) . we have 



so that 



dt dv 



TfC"(f)^- + fv 



(185) 



(186) 



, (187) 



where we have introduced the material derivative d/dt = 
d/dt + v ■ d/dr — W<P ■ d/dv in the advective term. This 
equation is more general than Eq. (|159[) . It shows that, 
for a given free energy, we can introduce an infinite class 
of NFP equations where £(r,v, t) is an arbitrary positive 
function of position, velocity and time. In particular, it 
can be a function of /(r, v, t). If we set £(r, v, t) = £g(f)/ f 
we recover Eq. (|159p . We can also write Eq. (|187|) in the 
alternative form 



df 

dt 



d_ 

dv 



£(r,v,t)[T 



1 



dv ~*~ C"(f) 



(188) 



where ^(r, v, t) is an arbitrary positive function of posi- 
tion, velocity and time. If we set £(r, v, t) = £h(f) we 
recover Eq. (| 1 59j) . These two alternative forms (|187| and 
(|188p were given in [33]. On the other hand, the general 
structure of Eq. (|185[) implies an iJ-theorem for the free 
energy (F < 0). The derivation is similar to that given in 
Sec. 12.10.11 



4.7 Equation of state 

The stationary solutions of the nonlinear Kramers equa- 
tion (|159p arc of the form / = /(e) with /'(e) < where 
e = i> 2 /2 + ^(r) is the energy of a particle. The func- 
tion / is determined by the convex function C accord- 
ing to Eq. (|176[) . Therefore, at equilibrium, the density 
p = J fdv and the pressure P = 2 I f v2 dv can be ex- 
pressed as p — p(<P(r)) and p = p(<P(r)). Eliminating 
the potential ^(r) between these expressions, we obtain a 
barotropic equation of state p = p{p) where the function 
p(p) is entirely determined by the convex function C(f). 
Furthermore, the condition that the distribution function 
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is a function / = /(e) of the energy alone implies the 
condition of hydrostatic balance. Indeed, we have 

Vp=-( f'(e)V<Pv 2 dv = [ ■ vdv 

a J d J dv 

= - W J fdv = -pV<P. (189) 

The condition of hydrostatic equilibrium can also be writ- 
ten p'(<P) = — p(#) or 



P'(p) 



(190) 



Let us introduce the free energy functional 

F[ P ] = Jpfyr d P >dv + \j P$ dr - ( 191 ) 

In Appendix [D] we show that this functional of p can be 
deduced from the free energy functional F[f] given by Eq. 
(|166[) by using the relation (|177[) valid at equilibrium. Fur- 
thermore, we show in Appendix lEl that the minimization 
problem (|181[) is equivalent to the minimization problem 



min {F[p} | M[p] = M}, 

p 



(192) 



where p(r) is the density profile corresponding to the dis- 
tribution function /(r, v). This equivalence considerably 
simplifies the study of the stability of a steady state of the 
NFP equation (|159[) . A critical point of F[p] at fixed mass 
satisfies the condition of hydrostatic balance. Indeed, writ- 
ing SF - a5M = 0, we have f p [p'{p')/p']dp' + <P - a = 
implying Vp + pV4> = 0. On the other hand, taking the 
second variations of F and using Eq. (|190[) , the condition 
of stability can be written 

S 2 F[Sp] =-\[J { j$) dr - J S P S<Pdr } * 0, (193) 
for all perturbations Sp that conserve mass. 



4.8 The strong friction limit 

In this section, we shall derive the generalized Smolu- 
chowski equation from the generalized Kramers equation 
in the strong friction limit £ — * +oo. The general case 
where both the diffusion coefficient and the friction coef- 
ficient depend on the distribution function is treated in 
[5l] by using a Chapman-Enskog expansion. Here, we re- 
strict ourselves to the generalized Kramers equation with 
constant friction coefficient 



dt 



df 
Or 



dv 



d_ 
dv 



TfC"(f)^- + fv 



with Eq. (|194| . Defining the density and the local velocity 

by 



(194) 



In that case, it is possible to develop a procedure simpler 
that the Chapman-Enskog expansion (see [71]). Let us de- 
rive the hierarchy of hydrodynamic equations associated 



P= fdv, pu = / fv dv 



(195) 



and integrating Eq. (|194[) on velocity, we get the continuity 
equation 



dp 
dt 



V • (pu) = 0. 



(196) 



Next, multiplying Eq. (|194j) by v and integrating on ve- 
locity, we obtain the momentum equation 



i(p u i) + ir(p u ^ u o) 



dt 



dxj 



dl% 
dxj 



d<P , 

- p^. — ipui 

oxi 



where we have defined the pressure tensor 
Pij = / fuiiWj dv, 



(197) 



(198) 



where w = v — u is the relative velocity. Using Eq. (|196l) , 
the momentum equation can be rewritten in the form 



dui 



dui 



P[ -dt +Uj dx, 



dxj 



d$ _ 
P dx~ ^ 



(199) 



We now consider the strong friction limit £ — > +oo 
with fixed T. Since the term in parenthesis in Eq. (|194[) 
must vanish to leading order, we find that the out of equi- 
librium distribution function /o( r j v ?^) is given by 



C'(fo) = -P 



+ A(r,t) 



+ o(r 1 ), 



(200) 



where A(r, t) is a constant of integration that is determined 
by the density according to 



p(r,t) = / f dv = p[X(r,t)}. 



(201) 



Note that the distribution function fo is isotropic so that 
the velocity u(r, t) = 0(£ -1 ) and the pressure tensor 



Pij = pdij + 0(£ ) where p is given by 



p(r,t) = i J f v 2 dv=p[X(r,t)}. 



(202) 



Eliminating A(r, t) between the two expressions pOlj) and 
(|202p , we find that the fluid is barotropic with an equation 
of state p = p(p) entirely determined by the function C(/). 
Of course, this is the same equation of state as the one 
obtained at equilibrium (see Sec. I4.7|) . Now, considering 
the momentum equation (|197p in the limit £ — > +oo, we 
find that 



pu 



(Vp+pw) + o(r 2 )- 



(203) 
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Inserting this relation in the continuity equation (|196|) . we 
obtain the generalized Smoluchowski equation [33|74|54| : 



dp 

at 



= v 



(Vp + /jV#) 



(204) 



The free energy associated to this equation is 



F[ P ] = J pf P ~JT d P' d r + \ J P**- 



(205) 



It can be deduced from the free energy (| 1 66[) by using Eq. 
(gOPj) to express F[f) as a functional F[p] = F[f Q ) of the 
density (see Appendix [Dj. A direct calculation leads to 
the if-theorem 



F 



(Vp + pV<P) 2 dr < 0. 



(206) 



The stationary solutions of the generalized Smoluchowski 
equation (|204p are critical points of free energy at fixed 
mass. They satisfy the condition of hydrostatic balance 



Vp + p\7<P = 0. 



(207) 



From Lyapunov's direct method, we conclude that a steady 
state of the generalized Smoluchowski equation 



linearly dynamically stable iff it is a (local) minimum of 
F at fixed mass M. This corresponds to the minimization 
problem (fTM]) . 

The condition of hydrostatic balance (|20T|) only holds 
at equilibrium. In the strong friction limit £ — > +oo, the 
out-of-equilibrium distribution is of the form /o = /o(e') 
with e' = v 2 /2 + A(r, t). Taking the gradient of Eq. p32l) 
and using a procedure similar to that followed in Eq. (|189|) 
with A(r,i) in place of ^(r), we obtain 



Vp = -pVA. 



Since p = p(p) and A = A(p), this can be rewritten 



(208) 



P'(p) 



(209) 



so that the out-of-equilibrium chemical potential A(r, t) is 
given by 



Hp) = 



p'(x) 



dx. 



(210) 



At equilibrium, comparing Eqs. (|208[) and (|207[) . we have 
A(r) = <P(r) + a/0 and Eq. |200jl leads to Eq. ([176]) . 



4.9 The damped Euler equations 



The generalized Smoluchowski equation (1204p can also be 



obtained formally from the damped EulcrlH equations [33] : 

(211) 



f + V>u) = 0) 



■^ + ( u - V)u = --Vp- V<2>-£u. (212) 

The Lyapunov functional associated with the damped Eu- 
ler equations is the generalized free energy 



F\p,u] 



V -^FT dp'dv + \ I p$dv + I p^di 



(213) 



It satisfies an ff -theorem of the form (see Appendix |G|) : 
F = — I £pu 2 dr < 0. (214) 



A steady state of the damped Euler equations (|21ip - (|212p 
satisfies the condition of hydrostatic equilibrium (|207[) . 
Furthermore, from Lyapunov's direct method, it is linearly 
dynamically stable iff it is a minimum of the free energy 
(|213|) at fixed mass. This corresponds to the minimization 
problem 



min {F[p, u] 

p.u 



M[p] = M}, 



(215) 



The friction coefficient in Eq. ()212|) measures the impor- 
tance of inertial effects. For £ = 0, we recover the usual 
barotropic Euler equations. Alternatively, if we consider 
the strong friction limit £ — > +oo, we can neglect the in- 
ertial term in Eq. (|212p and we get 



P 



(216) 



Substituting this relation in the continuity equation (121 1[) , 
we obtain the generalized Smoluchowski equation 



^ = V- 
dt 



i(Vp + /A7<£) 



(217) 



The physical justification of the damped Euler equations 
(|21ip - (|212p is not clear. They can be obtained from Eqs. 
P96p ~ (|197p if we close the hierarchy by invoking a local 
thermodynamic equilibrium (L.T.E.) condition |33j . How- 
ever, the rigorous justification of this L.T.E. condition is 



Hydrodynamical Euler equations involving a friction force 
— £(|u|)u have appeared in various contexts for different rea- 
sons. We may mention, for example, bottom-wall friction in 
2D turbulence [75], frictional force of air on a turbulent soap 
film [70], effective dynamical friction in the process of violent 
relaxation for collisionless stellar systems [28] and Epstein or 
Stokes friction laws for the dynamics of dust particles in the 
solar nebula 1771. 
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not established, so this approach remains heuristic. Nev- 
ertheless, hydrodynamic equations (hyperbolic models) of 
the form (|211| - (|212[) have been proposed in the context 
of chemotaxis to describe the organization of endothelial 
cells [78179151] . They lead to the formation of filaments 
that are interpreted as the beginning of a vasculature. 
These filaments arc not obtained in the Kcllcr-Scgcl model 
(parabolic model) which leads to point-wise blow up or 
round aggregates [52I42| . 



5 Explicit examples 

In this section, we give explicit examples showing the pas- 
sage from the generalized Kramers equation to the gener- 
alized Smoluchowski equation in the strong friction limit. 



5.1 Isothermal systems: Boltzmann entropy 



If we consider the Boltzmann entropy 

S B [f} = - J flnfdrdv, (218) 



5.2 Polytropes: Tsallis entropy 

If we consider the Tsallis g-entropy 
1 



S q [f] 



9-1 



(/« - f)drdv, 



we get the polytropic Kramers equation 



dt 



d_ 
dv 



The stationary state is the polytropic distribution 



(g-l)ff ' 
fi e 

q 



(225) 



(226) 



(227) 



we get the ordinary Kramers equation 



df_ = d_ 
dt dv 



(219) 



The stationary state is the isothermal (Maxwcll-Boltzmann) 
distribution 



where \i is determined by the conservation of mass. The 
index n of the polytrope is related to the parameter q by 
the relation 

n=l + ^. (228) 

Isothermal distribution functions are recovered in the limit 
q — ► 1 (i.e. n — ► +oo). We shall consider q > so that C 
is convex. We have to distinguish two cases, (i) For q > 1, 
i.e. n > d/2, the distribution has a compact support since 
/ is defined only for e < e m = pq/[\q — l|/3] (it vanishes 
at e = e m ). For e > e m , we set / = 0. For q — > +oo, i.e. 
n = d/2, / is the Heaviside function, (ii) For q < 1, the 
distribution is defined for all energies. For large velocities, 
it behaves like / ~ v -( d ~ 2n )_ Therefore, the density and 
the pressure are finite only for n < — 1, i.e. d/(d + 2) < 
q < 1. Therefore the range of allowed parameters are 



(220) 



where A is determined by the conservation of mass. The 
equation of state is the isothermal one 



d+2 



q>l, 

< q < l, 



n > - (case 1), 
2 v ' 



n < — 1 (case 2). 



(229) 



(230) 



P = pT- 



(221) 



The distribution function (|227p leads to the polytropic 
equation of state (see Appendix [Hj 



In the strong friction limit, wc obtain the ordinary Smolu- 
chowski equation 



dt 



P = Kp 1 , 7 = 
For n > d/2 the polytropic constant is 



1 

n 



(231) 



(222) 



K = 



The corresponding free energy is the Boltzmann free en- 
ergy 

F[p] = T j p\npdr+^J p$ dv, (223) 

and the stationary solution is the Boltzmann distribution 
p = A'er^ (224) 
where A 1 = (2tt / ' [3) d / 2 A. 



AS d 2- 



a_ i r(d/2)r(l-d/2 + n) 



r(i 



-l/r, 



and for n < — 1, we have 



K 



AS d 2* 



t.ridwri-n) 



r{d/2 



—l/r, 



(232) 



(233) 



where A — (/3\q— ll/q) 1 ^ 11 x '. In the strong friction limit, 
we get the polytropic Smoluchowski equation 



3e = v- 

at 



pV$) 



(234) 
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The generalized free energy is the Tsallis free energy 
K 



F[ P ] 



7' 



- J (p^- p ) dv + l - J p$ dv, (235) 



and the stationary solution is the polytropic distribution 

(236) 



P 



7-1 
A — — — 



Other useful relations valid for polytropic distributions are 
given in |80| . We note that a polytropic distribution with 
index q in phase space yields a polytropic distribution with 
index 7 = 1 + 2(q — l)/[2 + d(q — 1)] in physical space. 
In this sense, Tsallis distributions are stable laws since 
the functions /(e) and p(<£) have a similar structure. By 
comparing Eq. (|2%|) to Eq. (ggg) or Eq. with Eq. 

(|225|) to Eq. (|235p . we note that K plays the same role 
in physical space as the temperature T = 1/(3 in phase 
space. It is sometimes called a "polytropic temperature" . 
Quite generally, we define the local kinetic temperature 



T(r) by §T(r) 



or equivalently p{v) = p(r)T(r). It 



is proportional to the velocity dispersion. Since, at equi- 
librium, p = p{<P) and p = p{<P), we conclude that T(r) = 
T[<P(r)) is a function of the potential <P. This is true for any 
barotropic fluid at equilibrium. Now, for a polytropic dis- 
tribution, using Eq. ([23T]) . we have T(r) = Kp' 1 ' 1 . Then, 
using Eq. (f236|) we obtain T(r) = K\ - ^<£(r) so that 
T is a linear function of <P with a gradient (7 — l)/7 = 
l/(n + 1) = 2(q - l)/[(d + 2)q - d] [80I8T] . 



temperatures T — > +00. On the other hand, in the com- 
pletely degenerate limit, the distribution is a step function 
/ = rjoH(e — cf) where tp = — A In A is the Fermi energy. 
This is the case when A — > 0, valid at low temperatures 
T — > 0. This limiting distribution corresponds to a poly- 
trope with index n = d/2. The distribution in physical 
space, obtained by integrating the Fcrmi-Dirac statistics 
(|240p on the velocity, can be written 



where /„ is the Fermi integral 



(241) 



(242) 



The quantum equation of state for fermions is given in 
parametric form by 



VoS d 22 



At high temperatures (f — ► +00) we recover the clas- 
sical isothermal law p = pT and at low temperatures 
[t — > 0) we get a polytropic equation of state p = Kp 1 with 
7 = (d + 2)j2 (i.e. n = d/2) and K - ^(^) 2/<1 M- 
In the strong friction limit, we get the fermionic Smolu- 
chowski equation (|204|) where the equation of state is given 
by (|243[) . The fermionic Smoluchowski-Poisson system has 
been studied in [28166183] . 



5.3 Fermions: Fermi-Dirac entropy 

If we consider the Fermi-Dirac entropy 

S FD [f] = - VQ /(-ln^+fl-^ )ln(l- 
we get the fermionic Kramers equation 



— ) \drd~v, 
VoJ J 

(237) 



'1 
dt 



d_ 
9v 



(238) 



The corresponding equation with normal diffusion and 
nonlinear friction is 



( 1 
dt 



d_ 



D 



+ /3/(l-//r ?0 )v 



(239) 



The stationary state is the Fermi-Dirac distribution func- 
tion 

/ = ^ J ^7, (240) 



1 + Ae^ ' 

where A > is determined by the conservation of mass. 
The Fermi-Dirac distribution function (|240p satisfies the 
constraint / < r/o which is related to the Pauli exclusion 
principle in quantum mechanics. The isothermal distri- 
bution function (|220[) is recovered in the non-degenerate 
limit / <^ ijq. This is the case when A — > +00, valid at high 



6 Conclusion 

In this paper, we have studied a general class of non- 
linear mean field Fokker-Planck equations [33] associated 
with a formalism of effective generalized thermodynam- 
ics (E.G.T.). We have given several physical examples 
of application and we have shown that NFP equations 
can provide generalizations of the standard Kcllcr-Scgel 
model describing the chemotaxis of biological populations. 
The main properties of these NFP equations are valid 
for a large class of entropic functionals, encompassing the 
Boltzmann, the Fermi-Dirac, the Bosc-Einstein and the 
Tsallis statistics. Indeed, the rich mathematical structure 
of these equations is almost independent on the precise 
form of the entropy. These results should therefore evi- 
dence which properties in statistical mechanics are specific 
to the standard Boltzmann entropy and which properties 
are valid for a larger class of entropies. The distinguished 
feature of the Boltzmann entropy is that it can be ob- 
tained from a combinatorial analysis, assuming that all 
the microstates are equiprobablc. However, cquiprobabil- 
ity of the microstates is a strong postulate and it is not 
clear whether it has a universal scope. The universality 
of the Boltzmann entropy has been criticized long ago by 
Einstein [84] who argued that the statistics applicable on 
a system depends on its underlying dynamics (see discus- 
sion in [85]). For example, the Boltzmann distribution can 
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be obtained from a stochastic process describing a classical 
random walk where the kicks have uniform sizes. However, 
different distributions emerge when the stochastic process 
describes a biased random walk where the kicks depend on 
the region where the particle happens to be. We think that 
this is the case in many physical systems. This results in a 
very complex geometrical structure of phase space (frac- 
tal, multi- fractal,...) leading to non-Boltzmannian distri- 
butions at equilibrium. Indeed, in such circumstances, the 
microstates are not equiprobablc since the system prefers 
certain regions of phase space rather than others. It would 
be interesting to derive the corresponding generalized en- 
tropies (Tsallis, Abe, Kaniadakis,...) directly from a com- 
binatorial analysis. For example, Tsallis entropy could be 
the natural entropy on a fractal phase-space. The selec- 
tion of the entropy demands a complete specification of 
the microdynamics of the system in agreement with the 
statement given long ago by Einstein [84] . 



Since C is convex, we have 

[<?'(/) - C'(fo)] (f - f ) > 0, 



(250) 



so that F < 0. On the other hand, using a procedure 
similar to that described in Sec. 14. 3[ it is straightforward 
to prove that the steady states of Eq. (|244|) are given by 
Eq. (|176p and that a steady state of Eq. (|244|) is dynami- 
cally stable iff it is a (local) minimum of F at fixed mass. 
The first two hydrodynamic equations associated with Eq. 
(12441) are 



dp 
dt 



V ■ {pu) = 0, 



(251) 



dt 



d<P 

dxi 



-pu,, 

T 



(252) 



A The isotropic BGK operator 

Let us consider a simple kinetic equation where the gener- 
alized Fokker-Planck operator in Eq. (|159[) is replaced by 
a generalized isotropic BGK operator. This equation has 
been introduced in Appendix A of [54] and we provide 
here some complements. The kinetic equation is 



dt 



Or 



V<2>- 



dv 



f - fo 



where /o(r, v, t) is given by 
C"(/o) = -0 



y + A(r,t) 



(244) 



(245) 



In the limit r — > 0, we can repeat the same arguments as 
in Sec. 14.81 and we obtain the generalized Smoluchowski 
equation 



dp 
dt 



V • [r(Vp + pV$)] ■ 



(253) 



or 



This equation can also be obtained from a Chapman- 
Enskog expansion [54] . Rigorous mathematical results have 
been obtained recently in [8"6"j . 



B Connection between dynamical and 
thermodynamical stability for nonlinear mean 
field Fokker-Planck equations 



fo = F i/3 



v 

T 



A(r,t) 



(246) 



with F(x) = (C'^i-x). The function A(r,t) is deter- 
mined by the density by writing p = J f$dv = p(A). First, 
we show that Eq. (|244[) admits an H-theorem for the free 
energy (|166p . Recalling that the left hand side (Vlasov 
term) conserves the energy and the Casimirs, hence F, we 
find that 



F 



TC'(f) + — 



f-fo 



dvdw. 



(247) 



Using Eq. (|245[) and the fact that J fodv = J fdv, we 
have the identity 



TC'(f ) 

■0 



v 
~2 



f-fo 



dvdw 



A(r, t)— — —drdw 



Therefore, we can rewrite Eq. (|247[) in the form 



T / [C'(f) - C'(fo)} 



f-fo 



drdv. 



(248) 



(249) 



Let us consider a small perturbation Sp around a station- 
ary solution of the nonlinear mean field Fokker-Planck 
equation ([5]). We write the time dependence of the per- 
turbation as Sp ~ e xt . It can be shown that A is real, so 
that the perturbation is damped exponentially for A < 
(stable) or increases exponentially for A > (unstable). 
The linearized Fokker-Planck equation can be written 



A<5p: 



-V • S3. 



(254) 



If we multiply Eq. (|254p by C"(p)5p and integrate the 
resulting expression over the volume, we get 

A / C"{p){S P ) 2 dY= [ S3 ■[C"{p)V5p + C'"{p)5pVp\dv, 



(255) 

where we have used an integration by parts in the r.h.s. 
Now, the linear variation of the current (|10p around equi- 
librium can be written 

- °— = h'{p)8 P Vp + h{p)VSp + /3g'(p)5pV$ + (3g{p)VS$. 

(256) 
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Using h(p) = C"(p)g(p) and the relation 

C"(p)V P = -0WP, 
resulting from Eq. (f25|) , we obtain 



(257) 



Dg(p) 



C"'{p)8pVp + C"(p)WSp + (3W5<P. 



(258) 



Then, Eq. (|255[) can be rewritten 



A J C"{p){5pfdv + [3 J S3- V5$dr = - J JJ^ dr - 



(259) 



Now, multiplying Eq. (]254p by 5<P and integrating over the 
volume we find that 



A J 5p5$dr = J 5 J • V5$dr, 



(260) 



where we have used an integration by parts in the r.h.s. 
Inserting Eq. (]2"STj]) in Eq. (|2"59"|1 . we obtain 



(S3) 2 
Dg(p) 



dr. 



(261) 



On the other hand, the second variations of the rate of 
free energy dissipation (|24|) around equilibrium are 



S 2 F 



0*J) 2 
pDg(p) 



dr. 



(262) 



and they are clearly negative. Inserting Eqs. (]262p and 
([40]) in Eq. ()261|) . we finally obtain the relation 



2\5 2 F = 8 2 F < 0. 



(263) 



This relation shows that a steady state of the nonlinear 
mean field Fokker-Planck equation ([8]) is linearly dynam- 
ically stable (A < 0) iff it is a minimum (at least local) 
of the free energy at fixed mass (5 2 F > 0). Therefore, 
dynamical and generalized thermodynamical stability co- 
incide. 



C Stability of the homogeneous phase and 
critical point 

We consider a homogeneous stationary solution p(r) = p 
of the nonlinear mean field Fokker-Planck equation (J8j) . 
The corresponding potential is <P(r) = <P = Up where 
U = J u(x)<ix. The dynamical evolution of a small per- 
turbation around equilibrium is given by the linearized 
equation 

^=V-[Dh(p)VSp + X9( P )V6^] 

= Dg(p) [C"(p)ASp + (3A5<P] , (264) 



where we have used Eqs. (120)) and (|22|) to get the second 
line. The perturbations Sp and d^> can be decomposed in 



Fourier modes of the form /(r, t) 
leading to the dispertion relation 



J/(k)e ,;kr e A(k )*dk 



A(k) = -Dg(p)k 2 [C"{p)-(3v{k)], 



(265) 



where we have used 5& = (2ir) d u(k)8p (convolution) and 
set v{k) = -(2ir) d u{k). The syste m is stable (A < 0) if 



C"{p) - (iv{k) > 0, 



(266) 



for all k and unstable (to some wavenumbcrs) otherwise. 
For a potential satisfying v < 0, the homogeneous phase 
is always stable. Otherwise, there exists a critical point in 
the problem. The homogeneous phase is stable for 



T>T C 



V(k)y 



C"{p) 



(267) 



where v{k) max is the largest value of v(k) achieved for k = 
fc*. On the other hand, for T < T c the homogeneous phase 
is unstable to wavenumbers satisfying v(k) > TC"(p). The 
growth rate of the mode k is given by Eq. (|265[) . This sta- 
bility analysis has been explicitcd for particular potentials 
of the form @ in [40187] . It has also been generalized to 
potentials of the form ((3|) in [88] . 

The condition of generalized thermodynamical stabil- 
ity demands that p is a minimum of the free energy F[p], 
given by Eq. |36]) . at fixed mass. Using Eq. f3T]) . the sta- 
bility criterion PD]) can be rewritten 

S 2 F=-^iJ C"{p){5p) 2 dr + 13 J 8p5$dr\ > 0, 

(268) 

for all perturbations Sp that conserve mass. We need there- 
fore to investigate the eigenvalue equation 



C"{p)5p + {35$ = \5p. 



(269) 



The homogeneous phase is a minimum of free energy at 
fixed mass (stable) iff all the eigenvalues A are positive. If 
at least one eigenvalue is negative, the homogeneous phase 
is an unstable saddle point of free energy at fixed mass. 
Solving the eigenvalue Eq. (]269p in Fourier space, we get 

(270) 



C"{p)-/3v(k) = A. 



The spectrum of eiganvalues is continuous and Eq. (|270|) 
determines the eigenmodc k corresponding to the eigen- 
value A. Therefore, the homogeneous phase is a minimum 
of free energy at fixed mass (stable) iff the l.h.s. of Eq. 
(|270p is positive for all k. This returns the condition (]266p . 
i.e. T > T c . Alternatively, for T < T c , there exists modes 
k such that A < implying S 2 F < for these modes. 
In that case, the homogeneous phase is a maximum or a 
saddle point of free energy at fixed mass (unstable) . Using 
the Parseval theorem, the second variations of free energy 
(12681) can be written 



5 2 F 



2(3' 



(27r) d / [C"(p)-0v(k)]\S P \ 2 dk, (271) 
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leading directly to the preceding results. We also check on 
this explicit example (homogeneous state) that the condi- 
tions of dynamical and generalized thcrmodynamical sta- 
bility coincide. As shown in Appendix IB") this is also true 
for inhomogeneous equilibrium distributions. 

Finally, let us show that T c corresponds to a bifurca- 
tion point. The general steady state of the nonlinear mean 
field Fokker-Planck equation ((5|) is given by Eq. (|29j) . Close 
to the bifurcation point, the inhomogeneous equilibrium 
density profile can be written p(r) = p + e(r) where p is 
the homogeneous solution and e(r) <C p- Substituting this 
relation in Eq. (|29|) and expanding the equation to first 
order in e, we get 



with 



C"(p)e(r) = -13 



c'\)dr'. 



(272) 



In Fourier space, this relation becomes 

C"(p)e(k) = pv(k)e(k). (273) 

This equation has a non-zero solution e(k) ^ provided 
that there exists a mode k = fc* such that C"(p) = Pv(k^). 
This precisely corresponds to the condition T = T c . There- 
fore, stable stationary inhomogeneous solutions appear for 
T < T c precisely when the homogeneous phase becomes 
unstable. For the BMF model [43], the stability analysis 
can be carried out explicitly and we find that the phase 
transition at T = T c is second order. More generally, there 
may exist other systems where inhomogeneous solutions 
are stable for T > T c . In that case, they are in "compe- 
tition" with the homogeneous solution. One of these so- 
lutions is stable (global minimum of free energy) and the 
other metastable (local minimum of free energy). This is 
usually associated with a first order phase transition. 



D Passage from F[f] to F[p] 

We assume that the distribution function /(r, v, t) is given 
by an expression of the form 



C'(f) = -13 



y + A(r,<) 



(274) 



This expression appeared at several occasions in our anal- 
ysis (see Sees. |4~3l H3 |4~8] and Appendix |EJ. Since C is 
convex, the foregoing relation can be reversed to give 



f = F 



(275) 



where F{x) = (C") _1 (— x). Since the distribution function 
is isotropic, the local velocity vanishes: u = 0. The density 
and the pressure are then given by 



P = / rfv = PiM: 



P 



fv 2 dv = p[X\. (276) 



The first relation determines A(r, t) as a function of the 
density p(r, t). Substituting Eq. (|275|) in Eq. ([275]) . we find 
that 



1 



P 



pd/: 



1 



V = — d+2-M/?A), 

(3 2 



+ DC 



g{x) = 2 — S d F{x + t)t—dt, 



(278) 



^ r-t-oo 

h{x) = -2?S d F(x + t)t*dt, (279) 

where Sd is the surface of a unit sphere in d dimensions. 
Eliminating A between the foregoing expressions, we find 
that the fluid is barotropic, in the sense that p = p(p) 
where the equation of state is entirely specified by C(f). 
We can now express the free energy (|166j) as a functional 
of p by writing F[p] = F[f]. The energy (|162[) is simply 
given by 



(280) 



On the other hand, the entropy (|164[) can be written 

2^5^ S f r+°° 
S = —pdJT- J dr J o C W + ^ ( 281 ) 

Integrating by parts and using C'\F (x)\ = —x, we find 
that 

od/2c p r+oo 

S = — I dr I F'(t + 0\)(t + (3\)t d/2 dt. 



d[3 d / 2 



(282) 



Integrating by parts one more time and using Eqs. (|277p . 
flgTgp and (HZSJ), we finally obtain 



S = J pdv + 13 J Xpdr. 



(283) 



Collecting all the previous results, the free energy (|166|) 
becomes 

F[p] = - J p(\ + l^dv+ l -J p<Pdr. (284) 

Finally, using the relation h! (x) = —g(x) obtained from 
Eqs. (|278[) and (|279[) by a simple integration by parts, it 
is easy to check that Eq. (|277[) implies 



so that 



A'(P) = - 



P'(p) 



n'2 



dp . 



(285) 



(286) 



P J P 

Hence, the free energy can be written more explicitly as 



(277) F\P] = j pf yrdp'dr + \ J p$dr. (287) 
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E Equivalence between the stability criteria 
of the generalized Kramers and 
Smoluchowski equations 



Let us consider the minimization problem (|18ip . We want 
to determine the distribution /*(r, v) which minimizes the 
free energy F[f] at fixed mass M[f] = M. To solve this 
minimization problem, we can proceed in two steps. First 
step: we determine the distribution /i(r, v) which mini- 
mizes F[f] at fixed density profile p(r) = J fdv. This gives 
a distribution /i[p(r), v] depending on p(r) and v. Substi- 
tuting this distribution in the functional F[f], we obtain 
a functional F[p] = F[fi] of the density. Second step: we 
determine the density p*(r) which minimizes F[p] at fixed 
mass M[p] = M. Finally, we have /*(r,v) = /i [p* (r) , v] . 

Let us be more explicit. If we fix the density profile 
p(r), the potential energy W[p] is automatically deter- 
mined. Therefore, minimizing F[f] = E[f] — TS[f] at 
fixed density profile is equivalent to minimizing F[f] = 
K[f] — TS[f] at fixed density profile, where K[f] is the 
kinetic energy. The distribution /i(r,v) that extremizes 
F[f] with the constraint J f dv = p(r) satisfies the first 
order variations SF + J X(r)S(f fdv)dr = 0, where A(r) is 
a Lagrange multiplier. This leads to 



C'(fi) = -0 



r 2 

v 



A(r) 



(288) 



where A(r) is related to p(r) by writing p = J f\dv. Since 
S 2 F = -T5 2 S = \T J C"{h){5f) 2 dvdv > the distribu- 
tion fi is a minimum of F[f] at fixed density profile. Now, 
we remark that Eq. (|288j) has the form (|274|) so that the 
functional F[p] = F[f±] is explicitly given by 



p<Pdr. 



(289) 



Therefore, we conclude that /*(r, v) = /i[p*(r),v] is a 
minimum of F[f] at fixed mass iff p*(r) is a minimum of 
F[p] at fixed mass. Thus, the variational problems (|181[) 
and (|192p arc equivalent for global minimization. 

We shall now show that they are equivalent for local 
minimization. A critical point of (|181[) is a local minimum 
of F[f] at fixed mass iff inequality (|182p is satisfied for all 
perturbations Sf that conserve mass. A critical point of 
P92p is a local minimum of F[p] at fixed mass iff inequality 
(I193P is satisfied for all perturbations Sp that conserve 
mass. In order to make the connection between the second 
order variations (| 182[) and (|193[) . the idea is to project 
the perturbation Sf on a suitable space and write 5f = 
Sf\\ + Sf± where 5f± is the orthogonal perturbation (this 
is a relatively general method that has been applied in 
different contexts; see |36|89|90j ). We can always write 
the perturbation in the form 



Sf = 8fn + Sf ± = 



Sp 



-f(e) + 6f A 



(290) 



where 5f± = Sf — 6fu ensures that all the perturba- 
tions arc considered. By construction, we have J 5 fdv = 



f Sfndv = 5p so that J Sf±dv = 0. Therefore, 6f± is or- 
thogonal to Sf\\ in the sense that 



Sf\\Sfi 



f'(e) 



dv 



5f±dv = 0. 



(291) 



Then, we readily obtain 
(Sf? 



dv 



(Sf±) 2 
(Sf±) 2 



dv 



dv 



(5f\\f, 

(Sp) 2 
ff'(e)dv 



(292) 



Now, a critical point of (|18ip is of the form / = /(e) 
with e = v 2 /2 + <?(r) and /'(e) < 0. This implies that 
p = J f{e)dv = p{$) and p'{$) = J f'(e)dv. Therefore, 
Eq. (|292p can be rewritten 



(Sf) 2 
/'(e) 



dv 



(Sfi 



/'(e) 



-dv 



(Sp) 2 

pHpy 



(293) 



Combining Eqs. |[2S3]), (TI52]) and {Tgg]) we finally obtain 



S 2 F[Sf] 



(Sfi 



-drdv + 5 2 F[Sp}. (294) 



If 8 2 F[5p] > for all perturbations Sp that conserve mass, 
then S 2 F[Sf] > for all perturbations Sf that conserve 
mass. Alternatively, if there exists a perturbation Sp* such 
that 6 2 F[6p*] < 0, by taking Sf* in the form (|2"9"01) with 
Sp = Sp* and Sf± = 0, we get 8 2 F[Sf*] = S 2 F[8p*] < 0. 
We conclude that /(r,v) is a local minimum of F[f] at 
fixed M iff p(r) is a local minimum of F[p] at fixed M. 
Thus: (|18ip (|192|) for local and global minimizations. 

There are several consequences to this result (see also 
the more detailed discussion in [90]): 

(i) We have seen that / is a linearly dynamically stable 
steady state of the generalized Kramers equation (|159p iff 
it is a (local) minimum of F[f] at fixed mass M[f] = M. 
On the other hand, we have seen that p is a linearly dy- 
namically stable steady state of the generalized Smolu- 
chowski equation (|44|) iff it is a (local) minimum of F[p] at 
fixed mass M[p] = M. According to the above-mentioned 
result ([TSUI <^> we conclude that: /» (r, v) = f x [p* (r), v] 
is a linearly dynamically stable steady state of the gener- 
alized Kramers equation iff p* (r) is a linearly dynamically 
stable steady state of the generalized Smoluchowski equa- 
tion. 

(ii) It can be shown that a distribution function which 
minimizes a functional of the form F[f] = E[f] — TS[f) 
(where T is a constant and S[f] is given by Eq. (|164p 
where C(/) is an arbitrary convex function) at fixed mass 
M[f] = M is a nonlinear ly @ dynamically stable station- 
ary solution of the Vlasov equation (see, e.g., the case of 
stellar systems [H]). Therefore, the generalized mean field 
Kramers equation (|159|> . which precisely solves this min- 
imization problem, can be used as a numerical algorithm 



8 We implicitly consider here the formal nonlinear dynamical 
stability in the sense of Holm et al. [91] . 
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to construct nonlincarly dynamically stable stationary so- 
lutions of the Vlasov equation. Indeed, if a distribution is 
linearly dynamically stable with respect to the generalized 
Kramers equation then it is nonlincarly dynamically sta- 
ble with respect to the Vlasov equation (but the converse 
may be wrong, see below). 

(iii) It can be shown that a density profile is a nonlin- 
carly dynamically stable stationary solution of the Eulcr 
equation with a barotropic equation of state iff it mini- 
mizes a functional of the form F[p] (where F[p] is given 
by Eq. 107])) at fixed mass M[p] = M (see, e.g., the case 
of barotropic stars [92] )■ Therefore, the generalized mean 
field Smoluchowski equation (|44]h which precisely solves 
this minimization problem, can be used as a numerical al- 
gorithm to construct nonlincarly dynamically stable sta- 
tionary solutions of the Euler equation. Indeed, a den- 
sity profile is linearly dynamically stable with respect to 
the generalized Smoluchowski equation iff it is nonlinearly 
dynamically stable with respect to the barotropic Eulcr 
equation. 

(iv) According to the above mentioned results, a dis- 
tribution function / = /(e) with /'(e) < is a nonlin- 
early dynamically stable stationary solution of the Vlasov 
equation if the corresponding barotropic gas with equa- 
tion of state p = p(p) is a nonlincarly dynamically stable 
stationary solution of the Eulcr equation. In astrophysics, 
this corresponds to the nonlinear Antonov first law |92| . 
However, the minimization of F[f] = E[f] — TS[f] at 
fixed mass M[f] = M (problem with one constraint) is 
just a sufficient condition of nonlinear dynamical stability 
with respect to the Vlasov equation. Thus, this minimiza- 
tion problem does not allow to construct all the nonlin- 
early dynamically stable stationary solutions of the Vlasov 
equation. A larger class of nonlinearly dynamically stable 
stationary solutions is obtained by maximizing a Casimir 
functional of the form S[f] (where S[f] is given by Eq. 
fT64]) ) at fixed mass M[f] = M and energy E[f] = E 
(problem with two constraints) [92]. A numerical algo- 
rithm solving this maximization problem is proposed in 
[33 . In case of "ensemble inequivalence" , these solutions 
cannot be obtained by minimizing F[f] = E[f] — TS[f] at 
fixed mass M[f] = M. Therefore, the problem with two 
constraints ( "microcanonical" ) provides a refined condi- 
tion of nonlinear dynamical stability with respect to the 
problem with one constraint ( "canonical" ) . 



F Extension of the Eddington formula 



In Sec. I4.7[ we have seen that a distribution function / = 
/(e) determines a barotropic equation of state p = p(p) 
and we have explained how to obtain it (some explicit 
examples have been given in Sec. [5]). We shall now con- 
sider the inverse problem: find the distribution function 
/ = /(e) leading to the equation of state p(p). This prob- 
lem was first encountered in astrophysics and solved by 
Eddington [93j . In astrophysics, a distribution function of 
the form / = /(e) describes a particular class of spherical 
stellar systems that are stationary solutions of the Vlasov- 
Poisson system. To any such stellar system, we can asso- 



ciate a corresponding barotropic star with an equation of 
state p = p(p) which is a stationary solution of the Euler- 
Poisson system. The problem is to find the equation of 
state p = p{p) corresponding to the distribution function 
/ = /(e) and vice et versa. A similar problem arises in the 
context of nonlinear mean field Fokker-Planck equations. 
Thus, we can adapt many results of astrophysics to the 
present situation. For sake of generality, we shall consider 
the Eddington inverse problem in a space of arbitrary di- 
mension d. 

Knowing the equation of state p = p{p), we can obtain 
the equilibrium density p(<?) by integrating the condition 
of hydrostatic equilibrium (|190|) . The problem now is to 
determine /(e) from the knowledge of p(^). Let us rewrite 
the density in the form 



f(e)S d v d ~ i dv 



(295) 



We shall consider two cases: 

(i) We first assume that the distribution function has 
a compact support so that / = if e > e m . In that case, 
the range of integration in Eq. (|295j) is restricted to v < 
y/2{e m - W) so that 



Taking e = \ 
v we obtain 



f{e)S d v d ~ 1 dv. (296) 
<P as a variable of integration instead of 



p(*)= f(e)S d [2(e-<P)Y 



de. 



(297) 



It is convenient at this stage to define ip = e m — $ and 
x = e m — e. In terms of these variables, we get 



p^)=2^S d f f(x)(^-xf 
Jo 



2 dx. 



(298) 



In d = 3, taking the derivative of Eq. 



1 dp 
V8tt dip 



:dx. 



we find 



io VW- x 
This is an Abel integral whose solution is 

1 d [ x dp 



(299) 



fix) 



dip 



^/8tt 2 dx J dip yjx - ip 



1 



x d 2 p dip 



1 



o dip 2 ^Jx - x[) 



dp 

dw 



ip=0 



(300) 



This is the Eddington formula [S3]. In d = 1, Eq. l|29"g]) 

becomes 



1 



V2 



zdx. 



(301) 
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Comparing with the previous case, we immediately have G Derivation of the //-theorems 



V2ir 



ff ^ 1 d f ( M ^ 

V2-7T ax Jo y/x-tp 

x dp djj 1 " 

dip y/X — lp \JX 



Finally, in d = 2, Eq. (|297|) reduces to 



p(#) = 2tt / /(e)de, 
and we get the very simple result 

2ir d<P 



(302) 



(303) 



(304) 



(ii) We now consider the case where the distribution 
function takes strictly positive values for all energies so 
that e m — ► +00. Then Eq. (|295|) becomes 

r+°° d-2 
p{$) = / /(«■)$, [2(e -#)]"»-&. (305) 

In d = 3, taking the derivative of Eq. ()305|) we obtain 

1 * /(«) „, (306) 



V8tt d<£ 7$ V7~^$ de ' 

This is an Abel integral whose solution is 

1 d f +co dp d$ 



1 



V8n 2 



+°° d 2 p d$ 



d$ y/<P-e 



lim 



<lf! 



In d = 1, Eq. ([305]) becomes 



V2 



/(f) 



:de, 



. (307) 



(308) 



and we get 



V2tt 



V2irde 
+oc dp d<£ 



p(<2>) 



lim 



/0 



d^ - e <P->+oo ^JW~~t 



(309) 



For d = 2, Eq. (|304[) remains unchanged. For example, 
the distribution function associated with the Fermi-Dirac 
statistics in physical space 



p(<P) 



00 



1 + e^+ a : 



is 



a (3 



87rcosh 2 [i(/3e- 



(310) 



(d = 2). (311) 



Let us derive the iJ-theorem (ffiij) for the NFP equation 
(|8 |) -([2 |) . The time variations of the entropy ([2Tj) and of the 
energy (| 1 7[) associated with an external potential are 

S = - J C'{p)%dv, E=U ext ^dr. (312) 

On the other hand, using Eqs. (fP8")) and ^j, the energy 
associated with a binary potential of interaction can be 
written 

E = \J p(r,t)u(\r-r'\)p(r',t)drdr'. (313) 
Its time derivative is 

E=\ l^(r,t)u(\v-T'\)p(r',t)drdr' 

+ i I p(r, t)«(|r - r'|)^(r', t)dvdv' . (314) 

Interchanging the dummy variables r and r' and using Eq. 
([2]), we finally obtain 



(315) 



Therefore, the time variation of the free energy (|23|) is 
given by 

F = E -TS = J{$ + TC'(p))^dr. (316) 
Using Eq. @ and integrating by parts, we obtain 

F = f J • (V$ + TC"(p)Vp)dr. (317) 



Inserting Eqs. |[TU] ) . ([20 ]) and ([35]) in Eq. (|3T7|) . we finally 
obtain the i?-theorem ([24]). 

Let us derive the if -theorem (|25|) for the NFP equation 
([8|)-([3]). The time variation of the energy given by Eq. (jT9j) 
is 

iff d<p , 



■/t***/^*- (3i8) 



Integrating the first term by parts and using Eq. it 
can be rewritten 

E = y (zi<2> - fc 2 <i> - Ap) 2 dr + J ^dr. (319) 

Inserting Eq. ([5]) and proceeding as above, we obtain the 
ff-theorem (1251) . 

Consider now the generalized Smoluchowski equation 
(|44|) . The time variation of the free energy (|47|) is 

F= / [<?> + r^d,'dr+^|dr (320) 
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Inserting Eq. (|44|) in Eq. ()320|) and integrating by parts 
we obtain the ff-theorem (|48|). 

Let us derive the if -theorem (| 168[) for the NFP equa- 
tion (|159p -([2|). The time variations of the entropy (|164l) 
and of the energy (|162p are 



df 

C"(/)^drdv, E 



v 2 \ df 



(321) 



For D = £ = 0, the NFP equation ([159]) -([2|) reduces to 
the Vlasov equation. The Vlasov equation conserves the 
energy and the Casimirs. Indeed, using integrations by 
parts, we have 



E = 



$ -v 



or OT , 



J (V<£ • v v • V$)fdrdv = 0, (322) 



and 



S = - 



<9r 



drdv 



S ' (C(/)V) " £ ' 



gg(/) 
drdv = 0. (323) 



Therefore, coming back to the NFP equation (|159p -(f2")l. 
the only contribution to the time variation of the free en- 
ergy (|166p comes from the Fokker-Planck current. Using 
Eqs. (|32Tj) and IpO]) . we have 



F 



+ TC'(f) ) ^fdrdv 



TC"(f) 



TC'(f) 
df 



Of 
8J 
dv 

3 drdv. 



drdv 



(324) 



H Polytropic equation of state 

For n > d/2 (case 1), the polytropic DF can we written 

/ = A(e m - e) n - d/2 . (327) 
The density and the pressure can be expressed as 

p = AS d Q ($), p=^ASM$), (328) 



ith 



ry /2(e m -*) 

Qk= I I e™ - * - — 



2 \ n— d/2 



^"^v. (329) 



Setting x = w 2 /[2(e m — <£)], we obtain 

Q fe = 2 (fc+d-2)/2 (ero _ #) n+fc/2 

x f (l-x) n - d l 2 x (k+d - 2 V 2 dx. (330) 
Jo 

The integral can be expressed in terms of Gamma func- 
tions leading to 



Q fe = 2 (fe+d-2)/2 (£m _ (? y l +fc/2 

r{{d + k)/2)r{\- d/2 + n) 
X r(l + k/2 + n) ' 



(331) 



Then, the density and the pressure can be expressed in 
terms of the potential <P as 

Aq , ^ nnd/2 ^ r(d/2)r(l-d/2 + n) 

p = ASct{e m - 0) 2 1 — — — r , (332) 

1 (1 + n) 



p=^L(e m - (333) 
n + 1 f (1 + n) 



Inserting Eqs. (TTHT]) . (Ti"63|) and (fTeS)) in Eq. (|324|) and 

integrating by parts we obtain the ff-theorem (|168[) . 

Let us finally derive the ff-theorem (|214[) for the damped 
Eulcr equations (|211|) . (|212|) and J2]). The time variation 
of the free energy (|2 13[) is 



where we have used the identity r(n + 1) = nr(n) to sim- 
plify the second expression. Eliminating the potential <P 
between these equations, we obtain the polytropic equa- 
tion of state (j23l"|) with K given by Eq. (|232|) . 

For n < — 1 (case 2), the polytropic DF can we written 

/ = A(e m + e) n - d ' 2 . (334) 

The density and the pressure can be expressed as 



F = 



f> 



i^dp'dr + £ + — -^dr 
p J p 2 J at 

du 



2 
pu 



at 



dr. (325) 



Substituting Eqs. (|2TT|) and (|2"T2"|) in Eq. (|g2"5|) and inte- 
grating by parts, we get 



F= pu 



V 



Y ) -fu-(u-V)u 



dr. (326) 



Using (u • V)u = V(u 2 /2) — u x (V x u) we finally obtain 
the result ([2T4|) . 



p = ASM$), p = -AS d R 2 ($), (335) 



with 



2 \ n-d/2 



e m + $ + — 



v k+d - x dv. (336) 



Setting x = v 2 /[2(e m + <P)], we obtain 

Rk = 2 (fc+d-2)/2 (£m + ^n+k/2 
r+oo 

x / {l+x) n - d ' 2 x (k+d - 2 V 2 dx. (337) 
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The integral can be expressed in terms of Gamma func- 
tions leading to 



Rk = 2 ( k + d - 2 )/ 2 (e m + $) n + k / 2 
r((d + fc)/2)r(-fc/2-n) 
r(d/2-n) 



(338) 



Then, the density and the pressure can be expressed in 
terms of the potential <P as 

p = AS d (e m + ^) m 2 d/a - i r ^ ) 2 r _ ( ;" ) , (339) 



AS d 
n + 1 



(e m + <£) 



n+l 



2^i «H (340) 
r(d/2-n) ' 1 j 



where we have used the identity r(n + 1) = nr(n) to sim- 
plify the second expression. Eliminating the potential ^ 
between these equations, we obtain the polytropic equa- 
tion of state (j2"3"T]) with K given by Eq. ([2"3"3"jh 

Note, as a final remark, that spatially homogeneous 
polytropic distributions are obtained by taking <?(r) = 
in the above expressions. 
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